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SUMMARY 

This report is concerned with the problem of determining 
the influence of solar heating on the low frequency, dynamic stability of 
a gravity gradient satellite. The scope of the effort includes the develop- 
ment of a computer program that uses a quasi-static approach for 
describing a satellite’ s unsteady orbital motion. Tied to this objective 
IS a requirement to explain the source of the observed anomalous unstable 
behavior of the Naval Research Laboratory’ s gravity gradient satellite 
164 (NRL 164). 

The basic equations describing the rotational dynamics of NRL 164 
are obtamed in terms of kinematic, dynamic, orbital mechanic and 
thermal distortion equations., The resulting derivation leads to a set of 
fifteen first order, nonlinear differential equations of motion. The essential 
feature of the quasi- static approach is that the satellite’ s long, slender 
booms reach their thermal equilibrium position solely as a result of 
thermal bending and thermal lag without consideration of the dynamic effects 
of boom vibration. A computer program based on this formulation has 
been developed and used to examine the effects of absorptivity, sun angle, 
thermal lag and hinge stiffness on the stability of the satellite’ s motion. 

It is concluded that within the confines of relatively narrow 
stability criteria, the quasi-static model provides a valid means of pre- 
dicting the anomalous behavior of NRL 164. The occurrence of computed 
unstable responses closely resembling flight data tended to confirm 
that the source of the instability is related to thermal distortion and is 
particularly sensitive to sun angle. 

If future gravity gradient spacecraft stabilization systems are to be 
considered, a quasi- static mathematical model of sufficient generality 
should be developed for determining optimum stable configurations. 
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INSTABILITY OF A GRAVITY GRADIENT SATELLITE 
DUE TO THERMAL DISTORTION 


By Robert L. Goldman 
Martin Marietta Laboratories 
Martin Marietta Corporation 


INTRODUCTION 


The apparent simplicity offered by the concept of gravity gradient 
control of an earth orbiting satellite (References 1 to 3) has had a stimu- 
lating effect on the design of several spacecraft configurations. From an 
analytical point of view it seemed clear that three -axis passive stabilization 
of a spacecraft could be achieved through the clever use of tip weighted 
extendable booms. If stability were possible, a desirable earth-pointing 
equilibrium attitude of a satellite could be attained solely by a judicious 
arrangement of these booms (e. g. References 4 to 6). 

However, discovery of actual flight instabilities on several three- 
axis gravity gradient satellites made much of the earlier theoretical analysis 
suspect. Such instabilities were seen in the flight data collected during a 
series of gravity gradient experiments conducted by the Naval Research 
Laboratory (NRL) and examined in Reference 7.* Here it was observed that 
the unstable satellite response was primarily characterized by low frequency 
rigid body oscillations dominated by large yaw motions and, in some cases, 
yaw inversions. The satellite* s behavior was clearly related to the angle 
between the sun* s vector and the normal to the satellite* s orbital plane (sun 
angle), a condition that led directly to the conclusion that a satellite* s thermal 
distortion properties play an important role in the instability mechanism. 

If a correct analytical model of this instability mechanism could be 
developed, it might be possible to design gravity gradient satellites in such 
a way as to eliminate the detrimental effects of thermal distortion. As reported 
in Reference 8, a specific analytical attempt was made to determine the in- 
fluence that thermal distortion might have had on the attitude stability of the 
Naval Research Laboratory* s gravity gradient satellite 164 (NRL 164) shown in 
Figure 1. As described in Reference 7, NRL 164 was unstable in yaw while in 
full sunlight. 

The analytical model selected was the IMP Dynamics Computer Program, 
Reference 9. Although the program was developed to simulate the dynamics 
of IMP class satellites, it retains sufficient generality to be applicable to 
NRL 164. Internal and external effects due to gravity gradient forces, solar 
pressure, magnetic torques, eccentricity, thermal bending, boom vibration and 
a single -axis libration damper can be simulated by the program. 


’'This report is reproduced as Appendix E. 
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Unfortunately the planned 
analysis of NRL 164 with the IMP 
Dynamics Computer Program was 
unsuccessful. In retrospect the IMP 
Program was not at all suitable for 
the attempted investigation. In 
addition to a number of unanticipated 
program errors, the method proved 
to be much too cumbersome for the 
specific problem under investigation. 

Boom vibrations, for example, had 
an adyerse influence on the compu- 
tational procedure, since they led 
to excessively long computer runs. 

An alternative, quasi-static 
procedure, based on the dynamical 
equations described by Hooker in 
Reference 10, was suggested in 
Reference 8 as a means of simpli- 
fying the analytical model and the 
resulting computational procedures. 

In essence, this approach, which is 
similar to that used by Kanning in 
Reference 11, would eliminate en- 
tirely consideration of a boom' s Figure 1 - Satellite geometry, NRL 164. 
flexural modes of vibration. 

Assumptions of importance in the 
quasi-static approach are that: 

1. The boom reaches its thermal equilibrium position solely as 
a result of thermal bending and thermal lag without any con- 
sideration of boom vibration; 

2. The mertial and geometric properties of the satellite are 
altered as the tip weights at the ends of massless booms are 
displaced by boom deformation; 

3. The sun line and thermal properties of a boom determine the 
magnitude and direction of a boom* s tip deflection. 

The quasi-static approach eliminates the numerical analysis problem that 
arises from boom vibration and is justified by the relatively large ratio 
between the satellite* s boom vibration frequencies and its gravity induced 
rigid body frequencies. Boom bending is considered solely in terms of 
static thermal deformation without the local dynamic effects of boom inertia. 

In the present report the problem of ascertaining the influence of 
thermal deformation due to solar radiation on the low frequency, dynamic 
stability of NRL l64 is re-exammed. In this case the effort is directed 
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toward the development of a digital computer program that uses the 
quasi-static approach for describing the rotational dynamics of the 
satellite. The principal disturbance mechanism considered is assumed 
to be due to geometrical and subsequent inertial changes resulting from 
thermal distortion. The essential objective is to determine under what ^ 
conditions the resulting re-orientation process becomes unstable. 


SATELLITE CHARACTERISTICS 


The basic geometry of satellite 164 is illustrated in Figure 1, The 
three-axis, two-body gravity gradient stabilization system consisted of 
three extendable, interlocked nonperforated, SPAR BI-STEM booms 
arranged in a symmetric pattern about the plane of the pitch-yaw axes. 

The primary body was made up of the payload, mam boom and front 
lateral boom (fixed to the payload); the secondary body consisted solely of 
the lateral damper boom. Nominally, the lateral booms were located in 
the horizontal pitch-roll plane. The damper boom was connected to the 
primary/body through a smgle-axis hinge mechanism that constrained 
boom motion to a vertical plane. The hinge provided hysteresis damping 
torques and torsion wire spring restoring torques. 

By skewing the horizontal principal axis of the damper boom out 
of the orbital plane, all motions become strongly coupled. Under these 
conditions, three-axis damping of the entire satellite is achieved by the 
single-degree-of-freedom motion of the damper boom about its hinge 
(Reference 12). 

Satellite 164 was stable throughout its initial period of eclipsing , 
orbits (passage through the earth* s shadow) and unstable in yaw sometime 
after its first excursion into a fully sunlit orbit. In fact, during the entire 
first passage through eclipsing orbits and first entrance into full sunlight, 
all attitude errors were small. The perturbations were confined to approxi- 
mately one-cycle-per-orbit oscillations in pitch and l/ 2 cycle-per-orbit 
oscillations in yaw. 

As the satellite, m full sunlight, approached the 0^ sun angle 
position* the amplitude of the l/ 2 cycle-per-orbit oscillations in yaw un- 
expectedly increased. The satellite rapidly became unstable with several 
yaw inversions and large arriplitude oscillations m both pitch and roll. 

After passing through the 0^ sun angle position, the instability ceased and 
was followed by a period of stable operation. 


^ o 

In References 7 and 8 this is the 180 sun angle position. The 

difference lies in whether the sun vector is defined as either being toward 

or away from the sun. In the present report it is always toward the sun. 
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SIMULATION 


For NRL 164 the vector dynamical equations developed by Hooker 
in Reference 10 are used since their application eliminates many of the 
difficulties common to the Lagrangian approach (e.g. Reference 4). With 
respect to the basic equations, the final quasi-static attitude of the two- 
body satellite, including associated kinematics, dynamics, orbital 
mechanics, and thermal distortion, leads eventually to a set of 15 first- 
order, non-linear differential equations of motion. 

Reference Frames 


In the present derivation the three orthogonal frames illustrated 
m Figures 2 and 3 have been chosen. The basic inertial frame, 

[a] ^ 

IS fixed in space with its origin at the earth* s center. 

Here 

a^ IS a unit vector towards the perigee of the satellite* s orbit, 

■ ^ 

IS a unit vector parallel to the satellite' s orbital angular 
momentum vector, and 

^2 "" ^ 3 ^ ^1 

It has been assumed that the satellite* s orbital plane is perpendicu- 
lar to the a^ vector. The local inertial frame is simply obtained by a 

parallel translation of the [a ] frame from the earth* s center to the 
satellite* s center of mass. 

The second frame is the local vertical, or orbital reference 

frame, 

[E] = 

which has its origin fixed at the satellite* s center of mass and moves 
with it along its orbital path. 
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Here 


IS a unit vector in the direction of the radius vector, R, 
from the earth* s center to the satellite* s center of mass, 

IS a-unit vector parallel to a^, and 
^ ^ 2 ^ ^ 3 ' 

The first two reference frames are related through the dVrect'ion 
cosine matrix [B], leading to the orthogonal- transformation \ i 


where 


^1 


®11 

^12 


^2 


®21 

®22 

®23 

.^3. 


_®31 ' 

®32 

®33 . 

®11 


■®32 ^ 

-sin 

4^. B 23 

®12 

= 

®31 

cos 

4'.b^3 


L^3J 


( 1 ) 




21 ’22 ,33 

and = true anomaly. 

The third frame is a body frame, 

[T] = [e^ ^2^3] 

that has its origin fixed in the satellite and is used in conjunction with the 
inertial frame and local vertical frame to define the satellite' s motion. 
The nominal orientation of the undeformed NRL 164 satellite in terms of 

the unit vectors e^ , and e^ is depicted in Figure 3. 

The rotational position of this body frame with respect to the local 
vertical frame is described through the orthogonal matrix transformation 
[ S] such that 


6 


( 2 ) 


' 

®1 


Sii 

®12 

cn 

1 



c 

®2 

= 

®21 

^22 

^23, 


^2 


L^J 


.Si 

^2 

S 33 ^ 


.S. 


where S are direction cosines. These direction cosines are derived 
- ij - 

in Appendix A in terms of the usual Euler angles p,a. and y (designated 
as pitch, roll and yaw respectively). 

Combining Equations 1 and 2 results in the relationship 


{?} = [SI [B] [Zy = (Ai'^ir} (3) 


where [A] is now the direction cosine matrix relating the body frame 
to the inertial frame. 


Kinematic Equations 

The kinematic equations in the present derivation consist of the 
set of direction cosine rate equations that relate the elements of matrix 
[A] to the inertial angular velocity vector, cTq , where, m the body 
frame, ■ . . - 


CO 


= + 


“ 2®2 


Use 3 


(4) 


and are the body axis components of inertial rotation. For a system 
of moving axes (see Reference 13) the inertial derivatives of vector {e} 
are- given by i ^ . 


where 




-u 


3 




-CO 


1 


-CO 


2 


CO 

1 


0 


(5) 


( 6 ) 
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Since in inertial space {a } = 0 

the same inertial derivatives of vector {e} obtained from Equation (3) 
results in the expression 

{j’} = [A]^{r} (7) 

Equations (3), (5) and (7) yield the identity condition that 

-[ngl (Al^ 

which when transposed becomes 

[A] = [A] [S2g] 


When expanded this leads directly to the desired set of kinematic equations 
given by the six first order differential equations 


“^11 

“ ■ “^IS^Z 


'^u‘^ 3 ■^13“l 

Ai3 

= ^ll“2 ■ ^12“l 

An 

= ^22^3 ■ ■^23“2 

^12 

^2l“3 ‘^23“l 

^23 

“^21*^2 " ■^22“l 


(8) 


Only six kinematic equations are necessary since the orthogonality of 
the direction cosine matrix [A] requires that 


A 31 

- A A 
“^12 23 

■•^13^22 


-^32 

= Aj3A2j 

■ ■^n'^23 

(9) 

^33 

= -^11^22 

■ •^12'^21 
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Dynamic Eqgations 


The complete dynamic equations for NRL 164 are developed using 
the results of Reference 10. The vector dynamical equations with con- 
straint torques can be written in the form 




^^00 ^ ^ 10 ^ ' ^*01 ^ ^ 11 ^ * ^® 1 ^ ^0 ^ ^1 


(10) 


8r t^io • “o + ^11 ‘ (“o + ’^^1^1 = h' 


Here a tilde over an element is used to identify a dyadic, and a dot over 
a vector denotes its time derivative with respect to inertial space. 



$ 

ij 


the set of reduced inertia dyadics derived in Appendix B, 


6 = the angular rotation of the damper boom about its hinge 

axis, 

gj = the vector direction of the hinge axis (see Figure 3), and 

= the torque vector (gyroscopic torques, gravity gradient 
torques, hinge acceleration torques, interaction torques, 
etc. ) as defined in Appendix B. 


Written together in matrix notation Equation (10) becomes 


®00 

1 

- 

0 


g; . E ; 

O 

A 

1 

rH 

i-H 

- 

'b 

. - 


«i -_^i 


where aQQ is the dyadic 


'00 


= ^00.+ *10 +-®01 + *11 


uu 

i 


(12) 


a^j and a^^ are the vectors 


^01 ^^01 ^11^ ■ h 


(13) 
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A, such that 


^10 = h • ^^10 + ^ 11 ^ 


(1-4) 


and a^^ is the scalar 


hi ' Si ■ ^11 ’ ®1 


(15) 


Introducing the new dependent variable, 
6 = A 


( 16 ) 


permits the vector dynamical relationship of Equation (11) to be expanded 
into the scalar form 



(17) 


At this point Equation (17) is premultiplied by [Iqq] so that, along with 

Equation (16 the final complete set of dynamic equations is given by the 
five first-order differential equations 

-1 - , ' 








2 




— 

I 



00 




3 



A 





• 


6 = A 



(18) 


Orbital Mechanics 

The orbital path of the satellite is derived directly from the solution 
of Kepler' s equation based on a spherical earth model. This solution is 
carried out under the assumption that perigee occurs at t = 0 at which time the 
orbit path crosses the inertial a^^ axis. Referring to Figure 2 the complete 
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formulation of the decoupled Keplerian orbit is given by the first order 
differential equation 


4 ; = 




(1 


(1 + € COS l^) 


( 19 ) 


where 

= ( 

0 

a , = semi -major axis 

e = orbital eccentricity 

|x = gravitational constant 

ip = true anomaly 

The magnitude of the radius vector, R, is determined from 




= mean orbital rate 


|Rl 


a(l-€^) 

1 + € COSlp 


and it therefore follows that 


|R I cos ip 
|R[ sin ip 

Thermal Deformation 


R 

X 



( 20 ) 


( 21 ) 


Thermal deformation of each of the satellite' s three booms is pre- 
sumed to be due solely to a thermal bending moment induced by solar 
heating. The resulting deformation alters the position of each boom' s 
tip weight and, in accordance with the quasi- static approach, changes 
the satellites inertial properties. As the inertial properties vary, a new 
set of principal axes appear, and reorientation of the satellite occurs. 
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Deflection of a typical 
long, slender boom exposed 
to solar heating is illustrated 
m Figure 4. The undeformed 
boom lies along the unit 

vector axis, U, with the sun* s 

unit vector, cr , tilted so that 
it makes an angle ^ with the 
boom. The thermal bending 
moment induced by solar 
heating is considered to 
cause a tip deflection , cr 


i 

Deflected boom 

\ Ai 

T„ 



r 

Yj 

Sun vector 



A 




= y-j.^ 


( 22 ) 


Figure 4 - Boom deflection. 


•whose direction F , is normal to U in a plane containing the sun' s vector. 
The magnitude of the tip deflection, y^, in turn, is a function of the boom' s 

thermal properties as well as the position of the satellite with respect to 
the earth' s shadow. 

Let the unit sun vector, cr , be defined such that 
~5- - sin a + cos a ^ 2 , 


where a = sun angle. Since the direction of the unit tip deflection vector, 
r , IS assumed to always be away from the sun it can be derived from 
the vector triple product relationship 

^ _ (o^x !T) XU ^ 24 ) 

l(rxiT)xul 

Similarly the sun angle on the boom, is obtained from the equation 


cos C = U • O' 


(25) 


where, in view of the restrictions of Equation (24) 

0< ^<11 
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The tip deflection, y^, for each boom is computed from the first - 
order differential equation 

'll. yT’^T^T^T^TO ( 26 ) 

where t~= thermal lag and y™« is the static thermal tip deflection of the 
boom. ^ 

The static deflection, y,j,Q> is assumed to be the tip deflection that 

the boom would attain if thermal lag were ignored. The computation of 
y,j,Q is based in part on the analysis described in Reference 14. 

For a seamless, thin-walled cantilevered-boom of circular cross 
section it has been found that y.pQ caui be approximated by the equation 

2 

> - A f r A fcos^l 

y^0= ^4— J 

where f = boom length, ft , and is the thermal constant 



OoeDJg 

4hk 


(28) 


whose elements are 

= ^absorptivity 
D = boom diameter, in. 

e = linear thermal coefficient of expansion, 

in. / in. -°F 

h = thickness, in. 

2 

= solar radiation intensity, B tu / hr -ft 
k = thermal conductivity, Btu/hr-ft-°F 
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On the other hand if the satellite is in the earth' s shadow, solar 
radiation can be ignored and 


yTo = 


0 


( 29 ), 


The sun model used for finding out whether the satellite is either 
in full sunlight or in the-earth' s shadow is described in Appendix C.'"’ 

The three tip deflection vectors Y,pg, and Ypjj for the*T 

main boom, side boom and damper boom respectively are determined 
from the relationships of Equations (22) to (29). It should be noted that 
the imposition of- Equation (26) requires the solution of three first-order 
differential equations. 

Computed changes in the magnitudes and directions of the three tip 
deflections are used to rederive the reduced inertia dyadics $ 

(see Appendix B) and subsequently to alter the inertia matrix ^ 

[Iqq] in Equation (17). 


Computer Program 

The equations of motion describing the quasi-static response of 
NRL 164 to solar radiation are fully represented by Equations (8), (18), 

(19) and (26). The 15 first-order nonlinear differential equations of motion 
and associated supplementary equations have been programmed for digital 
solution on an IBM 360 computer. The resulting FORTRAN IV program, 
entitled GGSAT, is described in Appendix D.' GGSAT consists of a main 
program and a series of subroutines designed to solve the equations of 
motion. 

A solution is obtained using Hammings predictor -corrector method, 
with starting values determined by a Runge-Kutta procedure. A plotting 
routine, compatible with a CalComp 210/ 665 plotter, permits the user to 
generate plots of calculated variations of satellite yaw, pitch and roll as ' 
a function of time. - > ■ 

RESULTS AND DISCUSSION 


The essential aspects of the results presented here cover a 
selective investigation of the effect of boom thermal bending on the stability 
of the .quasi- static model of NRL 164 in a circular orbit (e = 0)'. Particular 
emphasis has been placed on the influence of absorptivity, sun angle, 
thermal lag and hinge stiffness on the satellites responses in yaw, pitch 
and roll. 
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The nominal physical and thermal properties of the satellite as 
well as its orbital characteristics are summarized in Table 1. The 
booms were considered to be long, slender beryllium copper tubes with 
highly reflective, silver-plated outer surfaces. For the sake of 
simplicity the mass of the booms was ignored. 


Stability Criteria 


The linear stability of NRL 164 in the absence of solar heating 
can be inferred from the roots of the linearized equations of motion given 
in Reference 5. The roots are either real or occur in complex pairs with 
each imaginary pair corresponding to one frequency of the satellite' s 


normal modes of oscilla- 
tion. Despite the strong 
coupling between the de- 
grees of freedom repre- 
senting yaw, pitch, roll 
and damper boom motion, 
each frequency still repre- 
sents a predominant modal 
motion in one of these four 
degrees of freedom. In 
Figure 5 the computed fre- 
quencies (obtained from the 
imaginary part of each root) 
along with the predominant 
modal motion at that fre- 
quency are plotted for a 
variation in damper spring 
constant. Here all other 
properties are consistent 
with Table 1. For 
Kj^ >. 0008 ft-lb/ rad the 

roots are all complex with 
negative real parts so that 
four stable frequencies 
exist. As Kjq becomes 



Figure 5 - Linearized frequency variation. 


large the frequencies 
approach those of a configu- 
ration with a locked damper 
hinge. At the design two of the roots are simply negative real. 


The 


complex roots all have negative real parts so that the three remaining fre- 
quencies are still stable. For Kjj < . 00064 ft-lb/ rad, however, one of the 


real roots becomes positive, and roll, pitch and yaw motions diverge 
exponentially. For reference purposes, this divergence boundary delineates 
the theoretical region within which linear stability of NRL 164 cannot be 
attained. 
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The nonlinear stability of NRL 164 exposed to solar radiation has 
been deduced through a time domain solution of the quasi-static 
equations of motion. Although the selected variations in parameters 
were quite broad, it was evident early in the study that nonlinearities 
in the simulation and economic limitations on computer usage required 
a narrow definition of stability. 

The stability criteria eventually adopted assumed that the satellite 
was unstable if, after the model was initially perturbed by a set of "prime 
conditions" in yaw, pitch and roll, a yaw inversion occurred within the 

TABLE 1 

PHYSICAL PROPERTIES AND ORBITAL PARAMETERS 
OF NRL GRAVITY GRADIENT SATELLITE 164 


Lateral and damper boom tip mass, slugs. . . 0. 118 

Main boom tip mass, slugs 0. 159 

Payload mass, slugs 8,800 

Lateral and damper boom length, ft 35. 0 

Mam boom length, ft 60. 0 

Distance from payload C.G. to hinge, ft 1,375 

Lateral and damper boom diameter, in 0.25 

Main boom diameter, in.T 0.50 2 

Lateral and damper boom wall thickness, in 0, 14x10 ^ 

Mam boom wall thickness, in 0.20x10 

2 

Payload moment of inertia, slug -ft 4.0 

Yaw rotation of lateral boom, deg -30, 0 

Yaw rotation of hinge axis, deg 120. 0 

Roll rotation of main boom, deg 0. 0 

Pitch rotation of main boom, deg 0.0 

Null position of damper boom, deg 0.0 

Damper sprmg constant, ft -lb/ rad 0.714x10 

Damper damping constant, ft-lb-sec/ rad 0, 395 

Eccentricity 0,0 

Semi-major axis, km 7302.4 

Earth radius, km 3 ' * “ ' 2 6378. 2 ^ 

Gravitational constant, km / sec 3.986x10^ 

Mean orbital rate, rad/ sec 0. 1012x10 

Initial true anomaly, deg 0.0 

Sun angle, deg 0, 0 

Absorptivity 0. 13 

Linear, thermal coefficient of expansion, in/ m- F 0. 104x10 

Thermal conductivity, Btu/hr-in- F 4.167 

Heat radiation of source, Btu/hr-in^ 3.065 

Thermal lag, mm 0. 0 

Initial roll angle, deg -10. 0 

Initial pitch angle, deg . . . -30. 0 

Initial yaw angle, deg 30. 0 
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Figure 6 - Nominal satellite response, no sun. 

first four orbits of the earth. These prime conditions (established 
somewhat arbitrarily by trial and error) were V= 30^, p = -3 0°, 
q; -10° and S' = ^ = a= 0. 

Although this definition provided a way of evaluating the effect on 
stability of a parametric change it was strictly limited by its singular 
specification of initial prime conditions and the number of orbits 
examined. Despite these limitations, the criteria still led to some 
particularly revealing observations. 

Basic Characteristics 


The motions of NRL 164 in the absence of sunlight and mitially per- 
turbed by the prime conditions are plotted in Figure 6, This figure 
serves as an instructive indicator of the basic nonlinear dynamic be- 
havior of the nominal configuration and should be used in gaging the changes 
m response that are brought about by a parameter variation. As expected 
m this case, the satellite is very stable. Although yaw motions are 
initially quite large they damp out rapidly and the system is close to an 
equilibrium attitude soon after the fourth orbit. 

Since the study is part icularly concerned with yaw instabilities , 
some insight into large amplitude yaw responses without sunlight is needed 
in order to evaluate basic system performance. Two such responses are 
plotted in Figure 7 for an initial yaw perturbation of 30 . The difference 
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between the two responses arises simply from the change in the damper 
spring constant. As can be seen, yaw damping and frequency decrease 
as the damper spring is stiffened above its nominal value (the drop in 
frequency is also apparent in Figure 5). The appearance of both pitch and 
roll coupling with a yaw perturbation is a basic characteristic of this type 
of gravity gradient configuration and is necessary in order to achieve 
three-axis damping of the entire satellite. 


Frequencies associated with the nonlinear response were, • as 
expected, amplitude dependent. For a small initial perturbation, however, 
the nonlinear influence on frequency was negligible^ and the quasi-static 
frequencies agreed closely with the linear analysis (see,Figure 5). 


Effect of Absorptivity 


The present thermal instability supposition rests on the assumption 
that thermally induced changes in the satellite* s inertial properties due to 
boom bending may cause a yaw instability. To examine this possibility, 
satellite responses were com- 
puted for variations in tip 
deflections brought about by a 
change in absorptivity. Fig- 
ure 8 shows the linear variation 
of the main boom* s static tip 
deflection with increasing 
absorptivity for normal sun 
incidence. 


Since the outer surfaces 
of the booms were assumed to 
be highly reflective, a low 
absorptivity, = 0. 13, was 

specified as the nominal value. 
On the other hand, as indicated 
in Reference 15, data accumu- 
lated by GSFC suggests that 
factors other than absorptivity 
tend to influence thermal bend- 
ing,and large deflections may 
occur in spite of a low value of 
absorptivity. To compensate 
for these factors, an effective 
absorptivity of 0. 5 was con- 
sidered to represent the high 
end of the deflection scale. 


o 


o 

Q. 



Figure 8 - 


Main boom tip deflection for 
normal sun incidence. 


19 




The computed influence of absorptivity on the attitude motion of 
NRLi 164 IS illustrated in Figure 9. For these responses, the sun is normal 
to the orbital plane (sun angle at 0°),and thermal lag is zero. The plots 
reveal the existence of a yaw instability that is particularly sensitive to 
thermal distortion. For Q!q = 0.3 and 0.4, the calculated yaw attitude, 
though large, seems to be stable. Conversely for CX q = 0.45 and 0.6, a 

yaw inversion occurs prior to the second orbit,and the configuration is 
considered to be unstable. The onset of an inversion appears to be closely 
dependent upon yaw reaching an angle of -90°. 

Effect of Sun Angle 

The investigation of absorptivity led directly to the observation that 
the satellite' s computed response and stability could be strongly influ- 
enced by a small change in sun angle. This was previously noted in the 
flight data of Reference 7 and now tends to be confirmed by the present 
quasi- static analysis. 

The sensitivity to a change in sun angle is typified by the computer 
plots shown in Figure 10, which depict the yaw, pitch and roll motions of 
NRL 164 for sun angles of 0°, 10°, 20° and 30°. Here thermal lag is zero, 
and the magnitude of boom thermal absorptivity is 0. 5, a value which for 
the 60 ft. main boom and normal sun incidence results in a computed static 
tip deflection of approximately 5 ft. (see Figure 8). 

In Figure 10(a), with the sun angle at 0°, a yaw inversion starts 
almost immediately and, as previously seen in Figure 9, the satellite 
response is unstable. In Figure 10(b), however, by just shifting the sun 
angle to 10 the inversion is suppressed and the motion is surprisingly ^ 
stable. In fact, in Figure 10(c), a further increase in the sun angle to 20 
still leads to a stable coijfiguration. Finally, in Figure 10(d), at a sun 
angle of 30°,occultations ' begin to occur, and the satellite is again unstable. 
The latter instability is probably involved with the deflection transient 
that arises as the satellite enters the earth' s shadow. 

Effect of Thermal Lag 

Figure 11 shows the effect on stability of a change in thermal lag 
with the sun angle at 0° and a g = 0. 5. Previously , with no thermal lag, 
(Figure 10(a)) this condition was unstable. For the time span shown in 
Figure 11, the inclusion of thermal lag appears to have a stabilizing 
tendency. 


Occultations are indicated by the downward facing steps in the 
horizontal line across the top of Figure 10(d), 
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Despite this observatioA it is not at all clear whether or not 
certain critical combinations of thermal lag, absorptivity, sun angle 
and damper spring constant may still result in an instability. In this 
respect it should be noted that in Figure 11, as thermal lag is increased 
from 10 min to 30 min , the magnitude of the second peak in yaw 
increases substantially. A further increase in thermal lag to 40 min , 
Figure 11(d), however, has the opposite effect, and the yaw amplitude 
begins to decrease. 

Effect of Hinge Stiffness 

The comparison of most of the computed results with flight data is 
in itself quite remarkable. As in flight, the controlling factor in 
establishing NRL 164' s stability boundary appears to be most closely 
allied with sun angle. Hinge stiffness, however, also plays an important 
role in the results. 

Some indication of the sensitivity of the computed response to a 
change in hin^ stiffness is shown in the plots of Figure 12. Here the 
sun angle is 0° and cXf. =■ 0. 4. For the nominal hinge stiffness 
(Kj^ = . 000714 ft -lb/ rad) the yaw response is stable. By just increasing 

the stiffness about 10%, so that Kq is equal to . 0008 ft -lb/ rad, the yaw 
response becomes unstable. 

Figures 13, 14 and 15 illustrate how the computed quasi-static 
stability boundaries for sun angles of 0°, 10° and 30° were influenced 
by variations in hinge stiffness (damper spring constant) Kjj, and 
absorptivity, a q. The illustrations essentially summarize the outputs 
from a large number of computer runs. 

In these figures, the divergence region arises solely from the con- 
sideration that,in the absence of sunlight, the system becomes unstable 
for Kj) < . 00064 ft-lb/ rad. Conversely the unstable regions that are 
shown in the figures encompass areas of yaw instabilities that were 
deduced by applying the nonlinear stability criteria. 

In Figures 13 ajid 14 the stability boundary is not only favorably in- 
fluenced by a low value of hmge stiffness, but the stable region increases 
in area as the sun angle changes from 0° to 10°, For a sun angle of 30°, 
Figure 15, the satellite passes through the earth' s shadow , and the 
stability picture becomes quite confusing. Now the unstable regions are 
isolated into two pockets that restrict the onset of an instability to only 
certain limited combinations of absorptivity and hinge stiffness. 

In obtaining Figures 13 to 15 a number of especially interesting 
cases were found in which the satellite responded in a yaw spin. A 
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Figure 13 - Stability boundary, a = 0°. 



Figure 14 - Stability boundary, a = 10°. 



Figure 15 - Stability boundary, a = 30°. 
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Figure 16 - Yaw spin, Kj^ = .001 ft-lb/rad, a = 10^, 
Oq = 0.425. 


typical occurrence is illustrated in Figure 16, in which the computed 
response in yaw rotates through 360° in about three orbits. Similar 
phenomena also were seen in the flight data of NRL 164 and conceivably 
were caused by parameter combinations of the type studied in the present 
investigation. 


CONCLUSIONS 


A nonlinear analytical model and a corresponding computer program 
have been developed to study the possible influence of solar heating on the 
anomalous low-frequency, orbital yaw instability of the Naval 
Research Laboratory' s gravity gradient satellite 164 in a circular orbit. 
The model' s formulation was based on a quasi-static approach in which 
the deflections of the satellite' s long, slender booms were determined in 
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terms of thermally induced bending without consideration of the dynamic 
effects of boom vibration. 

Within the confines of relatively narrow stability criteria, it has 
been found that, under the quasi- static model, NRL 164 not only becomes 
unstable but, in a number of cases, responses were computed that closely 
resembled flight data. 

From a review of the results, it has been discovered that; 

1. The onset of a yaw instability is particularly sensitive to a small 
change in sun angle, with the least favorable sun angle at 0°. 

2. In order to obtain an adverse effect from thermal distortion it 
was necessary to assume a larger thermal deflection than might be 
inferred from a boom' s nominal thermal properties. Some evidence is 
available, however, to justify this assumption. 

3. Occultation of the satellite produces isolated regions of 
instabilities believed to be related to the abrupt changes in boom deflections 
that occur as the satellite enters or leaves the earth' s shadow. 


4. In most cases, an increase in hinge stiffness tends to be de- 
stabilizing, while an increase in thermal lag seems to be stabilizing. 

By comparing the data collected during the computer study, it has 
been concluded that the quasi- static approach provides a valid means of 
predicting the anomalous behavior of NRL 164. In retrospect, the analysis 
probably could have been used to anticipate the instability seen in the 
flight data of NRL 164. 


RECOMMENDATIONS 

If three-axis gravity gradient stabilization of a satellite through the 
use of long, slender booms is to be considered in future spacecraft 
designs, adequate consideration must be given to the influence of' 
thermal distortion on system performance. Assuming that such satellites 
are still of practical consequence, it is important that the designer have 
available a means of selecting configurations that are free from the 
thermally induced instabilities. 

Since the present investigation has been limited to NRL 164 in a 
circular orbit.it is difficult to reach any definitive conclusion as to either 
an optimum configuration or ^the influence of eccentricity on stability. A 
general interrelationship between analysis and design selection is there- 
fore necessary before further development of three-axis stabilization 
systems can be adequately considered. In this respect, the formulation of 
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a quasi-static model of sufficient generality is recommended as the 
analytic tool for selecting and evaluating the nonlinear dynamics of a 
hinged multibody satellite. 

Closely related to this suggestion is the need to obtain a much 
j better physical understanding of the behavior of a boom exnosed to solar 
heating. A review is needed of existing GSFC experimental data ^with a 
view toward ascertaining the thermal bending process of a boom. This 
advice includes the possibility of additional experiments aimed at, determin- 
ing the effect of sun angle on the boorn and end conditions on actual thermal 
deflection. 
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APPENDIX A 
BODY ROTATIONS 
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An Euler rotational sequence of pitch, roll and yaw has been selected 
in order to define the orientation of the body frame, { ^} , relative to the 
local vertical frame, ,{E}- The three rotational transformations^ are: 


1. Pitch rotation Ej 



X ,T 


Here c = cos, s = sin 
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The total transformation is given by the matrix product 



c(3c\ + sccsPsY easy -spc\ + sq;c(3sY 

-cpsY + sasPcY cacY spsY + sacpcY 

casp -sa cacp 
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where are now the desired direction cosines that provide the 
relationship between the body frame and the local vertical frame 

In a similar manner, the angular velocity vector, w q, can be 
expressed in terms of the Euler angles and orbital motion by the 
matrix relationship 
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APPENDIX B 


INERTIA DYADICS AND TORQUE ELEMENTS 


Derivation of the elements of Equation (10) and (11) in the section on 
Dynamic Equations has been carried out in numerous references and is 
quite straightforward. 

The inertia dyadics are obtained from the geometrical and inertial 
properties of the satellite through the equations 


where 


1 


$ = 
00 

$0 + 

{L 

01 ■ ^01 ^ ■ ^01 

$ = 
10 

-m 

{=01 


$ = 
01 

-m 


•^01 ‘ -^iiSii 

$ = 
11 

m 


• L.JJ 1 - Ljj Ljj} 




m 


M, 


M, 


M, 


M- 


M. 


M 

L 


01 


= inertia dyadic of main body about its center of mass 

= reduced mass, MqM^/M 

= mass of main body, + M 2 + M^ 

= mass of lateral boom tip mass 

= mass of main boom tip mass 

= mass of payload 

= mass of damper boom tip mass 

= mass of satellite, Mq + M^ 

= vector from center of mass of main body to 
hinge point 
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vector from center of mass of damper boom to 
hinge point 


unit dyadic, + ®2 ^ ^3 


To illustrate the mechanics 
of deriving the elements of Equa-* 
tion (B 1 ), consider Figure B-1 , 
which shows the main body and 
hinge of NRL 164 with the booms 
deformed due to solar heating. 

The hinge vector, 
is given by the equation 


, Main Boom 


Lateral Boom 


Lo,= - 




C G 

“Mam Body 


Payload 
Damper Boom 


IVIyt 

^ (Secondary Body) 

where Figure B-1 Geometry of 

deformed satellite. 

Pj = vector from hinge to deflected lateral boom tip mass 


vector from hinge to deflected main boom tip mass 


= vector from hinge to payload center of mass 



= Pi 

+ Li„. 

1 

1 

01 




Ro 

= P, 

+ 

2 

2 

01 




^3 

= ^3 

'*■ ^01 


Then the inertia dyadic for the main body can be written 
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(B3) 


3 

io ' = '2 M. (IR.I^ 'l -5 r^)+ ij 

1= I 

where 

= the local inertia dyadic of the payload about its own 
center of mass 

The critical element in the quasi-static approach is the inclusion of 
the vector displacements Y^g and Y of each boom due to 

thermal distortion in the determination of the vectors P^, P 2 and L 
In terms of the nominal characteristics of NRL 164 (see Figure 3), the 
deformed position vectors are given by the equations 


Pj = i jCos <f)^ e j ^ sin 4> j e 2 + Y^g 

P2 = i 2®i^ 1 ’ 2 ^1 ^2® 2 ^TM 

^3 ^ 3 ® 3" ■ 

= i^sin92Cos 6 e -i ^cos *T* 2 Cos 6 e 2 - f ^ sin 6 e ^ - Y,^j^ 


The torque elements on the right of Equation (11) have 
been derived in Reference 10 and for the present case can be expanded 
into the following expressions: 

Eg + E* = 3G[ p .c + • p 1 - 6 [r,xiu' 8l'‘"0''<*00 + *li) 


• “0 


+ m L [ Wq X - 6 [wq X . g^ + g^ x 


^11- 


Wq] 


+ 6 m EQ2x[ci)Q X (^q^1-‘22)^ 

+ 5^mLQx[ (gjxL^^)] ^ , 


(B5) 
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{B 5 -Cont. ) 


+ Gm[ Lqj X (1 - 3 p p ) • L + Ljj X (1 - 3 p p ) . 


^ [^^01 ^ ^ 11 ^ ’ + «o ^11 ’ ^1 ®1 ’*^11 " “ol 

g 1 ' E = 3 Gg^ • [jr X . ^ ] - . [ Uq X • Oq ] 

- 6 gi • [ ^11 . gi + “0 ^11 • % + gi X • Wq ] 

- « ?1 • (?1 X *11 • «il + ”8i • [Ljix[a„«e„xL:gj)]' 

+ Gmgj^ • [ ^ (1 - 3pp) • Ijqj ] 

- S*' -Kd® 

where p = unit vector from the earth' s center toward the satellites 
composite center of mass 


Gx “o 


(!-«)■ 


(1 + € COS »j> ) 


Cjj = hinge damper damping constant 
Kjj = hinge damper spring constant 


in terms of the prescribed reference frames: 


p - cos ijj a^ + sinijja^ = ®i + ^23 ® 2 ^ ®33 ® 3 


= cos ^ 


g j = - sin <[>2003 e ^ + cos ^2 ‘*^3 ®2 ^ <|>2tOj^-cos ‘1*2 “2^ ® 3 
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APPENDIX C 


OCCULT ATION 


When the effect of solar heating on the bending of the booms xs included, 
consideration must be given to that portion of the satellite' s orbit during 
which the sun is occulted by the earth. 


In Figure C- 1 , occultation 
begins when the satellite in its 
orbital path reaches the occulta- 
tion point, Q2> Since the sun is 

assumed to be in the direction of 
the vector, the distances 

and 03^3 occultation 

point are given by 


^ 1^2 ~ 


(Cl) 


020^= QjQ 2 sinct= iRicos ijjsinQ! 


where is the true anomaly of 
the satellite and cosi|/sino£ < 0. 

Let the angle between the 
two radials002 and OQ3 be de- 
fined by the occultation angle v , 
where in general 



t 

I 


sin V = cos iji sin a 


(C 2 ) 


Figure C -1 Location of occultation point. 


Here OQ2 = IrI and OQ3 = R^, 

where R„ is the radius of the earth. Now since Q~ must be perpendicular 
E ^ 

to 02^3 occultation it follows that occultation occurs when 


or 



(C 3 ) 


(C 4 ) 
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Furthermore from consideration of the geometrical relationship between 
the sun vector, o' , and the occultation point, Q^, it is clear that the 
satellite must be in the sunlight as long as 
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COMPUTER PROGRAM GUIDE 

The accompanying information provides a brief description of the 
utilization of computer program GGSAT. The guide includes the required 
input data and a program listing as well as a short sample of the program 
output. ' 

The program is capable of solving for the satellite* s unsteady orbital 
motion under the influence of solar heating. It was written for an IBM/ 360 
operating system with plotting routines compatible with a CalComp 
210/665 plotting system. 


Program Input 

Actual data input is submitted as floating point information on three 
cards. Starting in column 1, data on each card is read in with a 4E16. 8 
format. The data must be supplied in the following order: 


CARD 1 


ALPHA, deg 
EPS, dimensionless 

ALO, dimensionless 


PSCALE, dimensionless 


Sun angle, a* 

Eccentricity, c , of 
orbital path. 

Absorptivity, 

0 . 0 < < 1 . 0 . 

Plot scale factor 
(usually 1. 0). 


CARD 2 ALFAE, deg 

BETAE, deg 
GAMAE, deg 
T OW , min 

CARD 3 V(15), ft-lb/ rad 


V(l6), ft-lb-sec/ rad 


Initial roll angle. 

Initial pitch angle. 

Initial yaw angle. 

Thermal lag, t. 

Damper spring constant, 
Kj^, (nominal value = 

. 000714). 

Damper damping constant, 
Cj^, (nominal value = .395). 
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RUNUM, dimensionless Number used to identify 

case being run. 

PLOTR, dimensionless Plot oi^^dinate selection 

for 8 j X 12" -Cal Comp 
plots abscissa rang- 

ing from 0 to 440 min at 
40 min per inch. 

' PLOTR = 1.0; ordinate 

ranges from -30° tO’+50° 
at 10° per inch. 

PLOTR = 2.0; ordinate 
ranges from -90° to 
+270° at 45° per inch. 

Program Listing 

The FORTRAN IV source deck for GGSAT consists of a main program 
and nineteen subroutines. Listings for the main p'rogram and eleven sub- 
routines are provided below. Of the remaining eight subroutines, subroutine 
HPCG is available from the IBM Scientific Subroutine Package, while sub- 
routines PLOTS, PLOT, FACTOR, AXIS, SYMBOL, NUMBER and LINE 
are part of the basic CalComp Software. 
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C “v( i ) =M 1 ,ilUE BjCK am MA.OAKHeK bOuK TIP MASo 
*. Vi(<: )=Kc,KAlN -BDCr Ilf fAiS 

C V( 3) =Mj .hAYLlAui fASS 

C V(A)=U rSIuE BliCM ANC .DAPPER flOuP LENOTn 

<- ■ «(b)=L^,PAU BoCf ..EACIP 

t vltl^LJ.blSTpNLE FrCP PAYLGAL C.&. TC HlAGt 

C . V(7 ) = I 11 ,11 .PAYLCAD POlL InEkTIA 
t. ' ■ V(81 = I(2,21 .fAYLLAD flTCf iPERTlA 
C '' V(V MI < 3 ,3 1 .f-A YLLAb YA« iNtPTlA 
C wI10) = PHniYAW RlTmTICc, LE SIDE BObP 

t V(il)=PHl2,YAW RLTaTICa lF HINGE AxM 

t ; 'Y112) = lhtTAl ,RuLl 'RClAflLN OF MAIN BLOW 
C V( 1 3 ) = THtTA2.,PlTk.H RLTATIOn EF PAIN bOuH 

t y( 1A)=ALPHA iERC.NoL POSITION OF DAPPER dOJP 

V. v(15)=RD,DAPFEr 5PrI.\G ClNS TANT , FT-Lo5 /R AD 1 An 

C YllGMtO.OAPPEK lArPINC C CNS 1 AN T , F T-LB-S EC /R AD 1 AA 

c . 


0 

c 

c 

c 

c 

L 

c 

0 

k. 

V. 

c 

t 

c 

c 

c 

L 

E 

k. 

E 

E 

E 

E 

E 

E 

E 

E 

E 

E 

C 

E 


PASS IN SLuC_S 
lENGIH in FEtT 
iNtRlIAS IN SLJG-FcEl*«2 
kCIATIjKS IN DtCREfcS 

alpha = aUN ANoLEt JEEREESI- ANGLE bETPEEN SLN AND NCBPAL Tl LRblTAL 
Plane vecicrs 

ALFAt = INITIAL rOlL AnC lE I DEGREES 1 
oElAt = initial PITCH ANCLE (DEGREES) 

CAPAt = INITIAL YAk ANgLE (DEGREES) 

a2 = SEMI-MaJLR axis (KH) CF CRoITAl PAlh 

tPs = eccentriliiv op lRbital path 

eRrDS = cARTP rAlIeS (KP ) 

p'Si .= initial true ANOHALY' (DEGREES) 

FMki ' gravitational ECNSTANT ( KM 3/S EE ** 2 ) 

EPtGA = PEAK CRBITAL RaTE (RaD/SEC) 

ISuN = 1 , CENSILEr ERtll IN SuN Sri A gCr ( ECl 1 PSc D ) 

= c, NlN-EELIPSEg LRblT (SoNSHlNt) 
mLU = ABSCkPlIVlTY 

iCe = LINEAR THErPAL CLEFFICIENT OF EXPAnS ION , IN/ ( IN-F I 

l 1 = SIDE BuCP CIAPETER, IKEhES 

g 2 = MAIN BUCP alAMtTER, INEHES 

LA =DaHPE(< bCEM DIAPIEIER, INEHES 

Hi = SIDE BUCP THICKNESS, INCHES 

H 2 = MAIN BuCP THlCKNcSs, INCHES 

HA =DAMPEP BCEM TH I E K NE 5 S , I NC HE S 

iC = THERPAL EDNDLCTmTY,BTU/(nR.-lN.-F) 

hRS = HEAT RaDIATIEN CF THE SOURCE, BTU/ ( HR- 1 N*« 2 ) 

PSEAlE. = PECl jCALE falter 

ID« = thermal lag time, input IN MINUTES^ 

RUNUP = RGN NuPtER FOR iOEN'T I E I CA T I CN 


external Out ,hCn 

LCMHlN/ Block 1/ AN, BN ,a2, CP EGA, EPS, NN, ISUN, ALPHA .ERRDS 
EOPMlN/BEDLK,!/T (23g ) ,R0U ( ,:3 C ) , P 1 TEH ( 230 ) , Y AR( 23 C) , I ISUN (230 ) » 
C rX(280 ),RY (230 ) ,R (230 ) 

ECMKlN/BLDY/V( IE) 

ECMPlN'/CLNsTS/PI .TkOPI , R AD 1 AN , FMU , DE GR EE ,S AR ,CAR ,ChI ,NDIP 
EDPMlN/X1N2/aLHAE,BETAc ,CAMAE,PS1 ,UMBC( 3) 

EDPPLN/HEATC/ATS ,ATM ,ATD ,ATC 1(3) ,ATC2(31 

LIPENSICN PRPT(E1,gEFFNC(9),AUX(U,13),Y(1S),DERY(15) 

UIMENSION 1BED(2) 

UATA IBCu/'CrBL' ,'T '/ 


r> r' o o o r/ 


bAiA Blank/ • '/ 

I.ATA IBRft/'RLN •/ 
lIhENSICN 1BCP(2) 

Ulf-ENSICN cSR(<:3C) 
tlHENSlCN hC (91 
i/ATA IbCK/'RtLL*/ 
jATA loCP/'PlTC *,'h •/ 

uATA IbCt/'YAW •/ 

i/ATA NO/ I n't.!! 5 .lit til 7. lid, 119, 120. 121/ 

PI=3 .1^,139265 

(.All plCis 

lWuPi=.:.«Pi 

kAUUN = PI/18C.O 

0EoRtE = l .C/RACIAN 

AZ=7;C^ .93 

EH(J = .39etliE46 

cRKDS=b378 . ItS 

l“BC(1)=0.0 

LHbC(2i =C.O 

LMbC (31=C.O 

PSl=L.O 

LiE L = C .0 

uELDlT=0.0 

9(1 )=0. 11759 

9(^1=0.1592 

9(31=8.7998 

9(91=35.0 

V(Bl=bO.C 

9(61=1.375 

v( 71=9.0 

V ( 6 1 =9 . 0 

9(91=9.0 

9( iOI=-3C.O 

9(111=120.0 

9(121=0.0 

9(131=0.0 

9(191=0.0 

TCE=.109fc-9 

Ul=.^5 

O2=.50 

09=.25 

r-l = l .9E-3 

h2=2.E-3 

C9=l .9E-3 

IC=9.167 

1-R5 = 3.065 


head sun angle alpha eccentricity tPS.ANL ABS0RPTI91 lY ,AL0 FkOH 
INPUT' CARDS AS" WELL AS P5CALE, PLOT SCALE FACTOR 
9(151,9(161 .RUNUM 

PLCTR=1. FCR Y RANGE CF -30 DEGREES Tu +50 CEoReES 

2. FOR Y RANGE- CF -90 OEGRtES TO +27C DEGREES 

1 READ(5,lC0G,£ND_=59q01ALPFA,EPS,AL0,PSCALt 
1000 FCRNAT(9Elb.Bl 

RE A 0(5,1 COO) AlFAE.BETAE .GAKAE.TCW 
READ(5,1C0019(151,9(16) .RUNUM.PLOTK 
ND1H=12 

1F(TUU.EU.0. IGU TO 2 
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IFITlH.M.IC .) Cl TC 55»CC 

N01M=lb 

lW=TLW*6C. 

LWi=i./Cw 

i. w1 = 1./FLlAT{ND1F) 

AR=ALPriA«RAClAM 
:>ak = >in(mR1 
lak=lc:> (ar ) 

ALFAbl=ALFAE 

DETAhI=BcTAE 

LAHAcI=GAM**E 

C lALClLATE uINEAR TnERMAL CuNjTANT FCR EACH BlCM 

ATC= (3.«ALa^TCc»hR5 ) /TC 
AT!i=ATC®Ll/Hi 
ATM = ATC«L2/H’£ 

ATD=ATC«LA/FA 

mTlI ll )=MTi« 

aTC 2(1 J=lATS«V(A))/3. 

MTCll2)=ATMo.5*Vl5)«w(b) 

«TL2 12 ) = (ATK^V(5 ))/3. 

ATl 1 (3 ) =ATU« .5«V (A ) (^1 ) 

ATl2i3) = IATC>>V(A))/3. 

L 

uP It l = i,<;3C 
1 ( 1 )=0.0 
KPLL(n=c.o 

j'lTChi n=o.c 

tAW( i)=0.0 
i^X ( I) =0 .0 
nVl I )=0.0 

K(i)= 0 .c 

lU CCMiNoE 
L 

LHfcGA=lFMU/(AZ*^3) 

FERIlD=(2.*P1/uPcGA) 

HRMTll)=C.C 
HRKT i2 )=^j6aCC. 
hRMT(3)=120.C 
pPHT lA ) = .1 
L 

CAlL UPoKDEPtAU 
Y(i)=DcPtNU(7) 

Y(^)=DtPcNu( c) 

Y(3)=0cPcNu(^) 

Y( A ) =DtLuCT<hAu IAN 
t(5 )=DtL»RACiAN 
Ylo)=DcPcK0( i) 

Y(7)=DtPtN0(^) 

Y(b)=DtPcNU(3) ' 

Y(^>=0kPcNu( A) 

YUC)=uEFEisD (5) 

Y(i1)=UEPEaC(6» ' - ^ 

YU2)=PSiOKA^UK ' ’ 

Y(i3)Bu«0 

Y( aA^=o.C ^ 

1 ( i 5 i *U »L 

V 

uO laS J»1 ,NlIM 
I5i lEkV(J)>.l 
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lRUN = ftUNl.H 
«iRiTt(6f<i999IIKUN 
^^99 fOHHAT(lhl,6CX, 'RUN'.IA,//) 

kRlTt(6,50Q0)(V(w) , lfDl,02,hl,H2 
hRlTt(6fS0iC)(W( J) 

hRlTt(6>bl00 J £P:>,AZ>EKR1;S,PSI . FHU.CREGA .ALPHA 
hRlTt:(6t5020)AL0.TLE»TC.|-R:>.T0W 
hRlTt(bi3150 ) ALFAEl .BETAEI .OAHAEI 
KRlTt(6.51eC) 

C 

AN=0.0 
bN =1 .0 
ftN = l 

riRlTE(6.5200) 

L 

call HPCC{PRPT,t,OtRlr,hCIH, 1 hLF.FCN.lUT .AUX) 

C 

AN=NA-1 

NNi=NN+l 

NN2=AN42 

T(AN1)=0.0 

1(AN2)=4C.0 

SS=A5. 

lL = 3. 

FV=-jO. 

LV=1C. 

IFIPLDTR.EU.U )Cl TD 136 

5S=2t5. 

lL = 2. 

FV = -‘.0. 
lV=A5. 

15b KGlL(Nn1 ) = FV 
kGlL(NN 2)=UV 
HlTCh(rtM)=FV 
HTCHNN2)=C^ 

YA(J (NNi )=FV 
YAh(NN2>=DY 
ESP(NN1)=FW 
tSP(ANt)=Dv 

L 

UO 160 1=1. NA 

iF(( YAk>(n .tC.PLAAK ).AAl. (PlTCHd ) .EO.BlANK))GO TU 157 
nRlTc(o.53LC )T(1 ) .H ( i) .RX( 1 ) .RY( 1) ,RLLL( 1) .PITCH! I ) .YAt. ( I) . 
C nSLMl) 

lG TU 16C 

157 nRlTt (b .3301 )T1 n.A ( 1) .RM n .RY ( 1) .RULl! 1) . IlSUN! I) 

5;01 rCKHAT(6X.AFl0.2.2PF12.2 .2AX.I5) 

Ibu CGNTlNuE 

CAlL PLn(l7..-15. .-3) 

CALL PLOTIO. .1.5.-3) 
iFlPiCALt.EC.l.OlGu 1C 165 
call FaC TCA (FSCAlE ) 

ibi lAlL AXI6(C. .0. .13HT1M6 . M I NOTE 5 . - 1 3 . 1 1 . .0 . , C . , A C . ) 
call Aa I3(u . .0 . . IAHAAGlE . 0E0RLE3 . 1 A.6 . .90 . .FV ,CV ) 

CALL PlCTIo. .ZL. 3) 
lAlL PLCmi ..ZL.2) 
iPACc=0. 
iFPA=0 

yf 2yO 1=2.9 

iPMCL = 6PACE + PEKlLD/ (12 l .«AU) 


I 


D-6 




call Pt.C1(SP«Ct,C.i3) 

YLEN = .i25^rtCCn 
CALL PLCKiPACfc.ILkN’.Z) 
lF(HLD( I fit .tO.Ct OC TL 190 
lFPA'*IFPA + i 
YP»SPACE-.125 

call SYHeDLl YP, K.AC(1FPM.0.,-1) 
cO TL 200 

19U 1F(I.NE.2)(>C TO 191 
XPACE»l>PACl; 

*C = Xt-AoE + .2 

CAlL SYMoOL(XPaCE,.8C,.1C.1.C.i-1) 

CAlL SYPBCcl>0,.75..i0,IbCK,0.,<i} , 

call SYHbOL(XPACE. .55, .1 Ci2>0.i-1) 

CALL $YHBDlIxC,.50, .10,IECF,0.,8) 
call SYHbCL(XPACE,.3C,.lC,3,C.,-i) 
call SYHdOl(xO, .< 5, .lC,ItCY,C.,AI 
xPC»xPAGk-.75 

LAuL SYnBQLlXPu,b.,.l,lBKN,0.,<>> 

XX=XHG».9 

CALL NoMbERtXX.t .,.1,RCALH,0.,-1) 

191 IFU.KE.h) GC TC 2CC 
aPAGe = SPACE- .25 

CALL SYHB0UxPAGE,.2S..lA,lBCD,0.,b> 
bOO CONTINUE 

CALL L1NE(T ( 1) ,RLLL( 1) ,Nh ,l.-^b,l) 
call LlNeCT ( 1) .PITCHU ) ,NN ,1,*8,2) 

CALL LlNE(T(i),YAM(ll.NK .l,-»8,3) 

OD 3u K>1,NN 

cSP(K)= iS + 5.*bLLAT(llSCMK)-l) 

3u CONTINUE 

CAlL L1NE( nn.ESPd ),NN,1,4C,11 ) 

CO Tu 1 

2u FOMMATI//, (AX,JF17.7J) 

21 FCaHAT ( //, (AX, aFl7 .7 ) ) 

2b F0KMttT(//,(AX,AFl7.71) 

2j E0RHAT(/y,iAX,fl7..7J ) 

SCOU F0kNAT(/.20X,>Ml ,S10E bCCH AND HA.OAHPER BOON TIP M ASS ° , F 10 .5 , 1 X , 1 
C* SLCGS'/ 2 

CbOx, ‘Pb.FAlN BbCF TIP PASS F9 .A_i2X SLCGS ' / 

CbOX , 'H j.PAYLLAu pass =',F9.A,2X,' SLUGS'/ 

CbOX, 'Ll.SIUE BUCP ANL lA, DAMPER BOOM LENGTH =',F7.2,AX,' FEET*/ 
CbOx , 'LbiMAlN BuCM LENGTH = * , F7 .2,AX, ' .FEET ' / 

CbOX, 'La.CIbTANCE FROM PAYLUAD C.G. TU HINGE =',F8.3,3X,' FEET'/ 
CbOx, '01 ,SIUE BuCP ANL LA, DAMPER BOOM 0 1 AMET EB= * , F7 . 2 ,AX , ' INCHES'/ 
CbO*. 'Db,MAlN BOOM UIAHETER =',F7^,AX4' INCHES'/ 

CbOX, 'H1,M0E,A.\D HA, CAMPER, BDUH THICKNESS =',1PE11.2,' INCHES*/ 
CbOX , 'HA,MAIN BUCM IH1C^^E5S = * , IPEl 1 .2 INCHES')' 

501U HORMATl/, 

C2CX , ' I { 1 ,11 .PAYICAC BOLL INERTIA =*,F7.2,AX,« SLCG-FT2 

C •/ 

CbOX, '^2,^) ,PAYU.Aa FITCH INERTIA ? • ,f , Alt , • SLUG-ET2 

C */ 

CbOX, '1(3,31, PAYLLAO YAk INERTIA =',F7.2,AX,' SLUG-FT2 

C •/ 

CbOX, 'P hIUYAr RCTATICN CF side BOOM =*,F7.2,AX,* DEG'/ 

CbOX, 'PHI2, YAk RCTATILN CF HINGE AXIS =*,F7.2,AX,' UEG'/ 

_CbOXi'THETA.I,KOlJ^ ROTATION OF MAIN BOOK =*jF7.2*AX.* UEG'/. 

CbOX, 'THETA2, PITCH RCIATILN CF MAIN BOOM =*,F7.2,AX,' OEG'/ 

CbCX, 'DELTA 2ERU,NULL PlSITIDN OF DAMPER BOOM =',F7.2,AX,' DEC'/ 
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C*^)a, 'K u,CAr.PtR SPRING LCi\STAUT 
Cu’/ , 

CiiOX, 'Cu,uA«PtP OAMPUG ClNSTANT 
CiEC/rAu'//> 

5i0u hOKMAT(/,2CX,'ECCeMiJICITY 
DiOXf 'StMl-HA.CR AXIS 
DtOX , 'EARTH RADlLS 
CaOX.'IMTIAL TkLe AACMALT 
DaCX , 'GRAkITATIlAAL CCAsTAA'T 
C ‘/iOx, 'OREGA ,HE^N LRtlTAL KATE 
DC’/ 

DiOx , 'Sum ANGLE 
C/) 

5C2U FORPmT(/, 

FtOX , 'AbSCRPT IVlTt 

FiOx, 'LINEAR THtRFAL CCEFUCIENT LF EXPANSICA 


= 't^PtU.x,' FT-LB/kA 

= • ,CPF8.3,3X, • FT-LE- 

=',F11.6/ 1 

= • fAPEH.i , ' KK ’/ i 

= 'iAPEli.C, ' KK’,/ 

= • ,APE11.2 , • LEG' ,/ 

= ',APEli.2, ' KK3/SEC23 
='.APE11.2,' KAl/SEA 

= • .LPF7.2,AX, • lEG',/ 


-*tF 7 * 2 tAXf/» 

= ’ iAPEll.2, ’ IN/IN-F ' 


F,/ 

F^OX, ‘Thermal CuACUCTiVlTT =‘,APEU.2,' uTU/hR-I 

FL-F •/ 

FtOX.'HEAT RALIaTION LF TFE SLUKCE = ' , APE 1 i . ’ BlL/nP-I 

FN2' ,/,20x, ‘T fErPAL LAG TlMt * ' , LP F7 . 2 ,Ax , ' 

GPIN • ,// ) 

5i50 FCkHAT(/,2UX,'iMTlAL RClL ANGLE =',FU.6,' D 

EtG* ,/ 

EtOX, 'INITIAL PITCH ANGLE DEG',/ 

E20X,'ImITIAL yaw angle =‘,F11.6,' OEu',/1 

5160 FOrMaTIIHII 

5<lOl; FCRWaT(//, 9X, 'j IPE ' , 8X , ‘A ' , 8 X , *RX ' , 8 X , ■ fi V • , 7X , ‘ROLL ’ 1 7X , *P 1 TCH • , 
CbX, ‘YA« ' ,5X , 'ECLlPiEC ,/ ) 

5 30u ECkPAT(AX,AF10.2 .2P3E12 .2i 15 I 
5500 WRlTt(s.560C)TQW 

5600 FDRHATdhl,' THERMAL LAG TIMt =',E16.8,' AnD IS LOT CF RANGE.' ) 
CD Tl 1 

5A00 CALL PLCr(12.,0.,9V9» 

STOP 

END 
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i^UfifUUT IhE LUTOt^ •lHLF,NClHtP(<HT) 

L 

cCHHLN/BLCCKi/ANtBN.^Z.ChEGA.EPS.NN.lSUN.ALPI-A.EKRDS 
COMHuN/BlOCKc/T ( 230 > .ROLL ( 2 30 ) . P ITCH (230 ) , Y AK( 230 1 , 1 1SUA(230». 
C kY(2301.KY(2^0).K(230) 

COHHLN/_CCNiTS/P I .T>iCPI .RAOIAN.FPO.OECREE.SAR.CAR 
LIHEASICA PRPT(5).Y(1>.CLRV(1).A(3,3>.BI3.3).S(3,3) 
lATA ALPlLO ,LErCL0,GAHCLC/i«0.0/ 

OATA Blank/ • . •/ 

4=X-AN*PhMT (2) 

1F( ZI5oC.20.lC 
10 PRMT(51=i.C 
OC Tu 500 
C 

2u Al=AA«>ll.-EPi**21 
LY12=CU5IY( li) ) 
iY12=SlN(Y( li) 1 
A2=(1.4EPS*CY12) 

H(hM = Al/A^ 

fcX(KN)=R(H«n«CYl^ 
kY(NM=R(NN)*SY 12 
l(NN)=X/t0.C 
1I5UMNN ) = ISLN 
o 

CALL EULtR( Y.A.e.Sl 

1F(S13, n.kO. 0.0 . And. 5(3. 3) .bO. 0.0)00 TO 50 
ti = -S(3.2) 

ALFAt=AR5lN(k) 

iFlAo5 ( AlFaF *D£CHE£ ) .CT . £9 . 0 .A«D .Ab 5 (AlF AE •OEGHE £ ) . LT . 9 1 ,0) 
loO To JO 
5 N=ClS(ALFaE ) 
iN'l = l./SN 
TE5Ti=J(3.j)«S.«l 
a = 5(3.1)*SM 
C 

i>ETAc=AR5IN(».) 

lF(TtSTl.LT.C.O)tE)At=-AKSlN(W)*PI 

C 

7E5T2 = 5 (2.2 )«SN1 
•> = A( i.2)«SM 
lAHAc’ARJIN ( k) 

1F( Tt5T2.LT.C.0)0APA£ = -A).5JN 0>)+Pl 
C 

MLFA£ = ALFAt«uE(>RkE 
tBTAt=bETAt<'LEGBtE 
OAHA£=oAMAE*LEoREE 
C 

aLPCLC=AlFaE 
bETOLD=B£TAE 
CAhCLD = CAliAE 
C 

(.0 T c 1 1 0 
C 

30 ALFA£=ALFA£ «CEoR £E 
ALPCLD=AtFAf 
lC TL lot 

50 lFiSl3.2)-l.C)bC.BC.6C 
60 iFlS(3.2)'»l.C)i0t..7C.U0 
7o ALFAt=9C.C 
«LPOLD=AlFaE 


I 


. CO TC IOC 
80 ALFAt=“9C.O 
ALPDLD=ALFAE 
100 KOLL(Kft)=ALPUO 
PITCF(NN)=BLaNK 
YAW(NM-8LAKK 
CD Tb 111 

110 CONTINUE 

ROLL (NN)=ALPLLD 

PITCh«NN)=bETOLC 

YAn(NN)=GAMOlD 

111 NN=NN+1 
AN = AA-^bN 

500 RETURN 
ENO 
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(,««»o»»o»»»o«o»»»»«»aSl.RKrLTlKt FCh»**** 

Subi^DuTiNt FC.«( A.V .OEkY) 

L 

l.r«»i,N/Bwni.Kl/AF .Cl-FOA .FHS.Nft, JSON , ALPnA , EkRUS 

tPhPoN/PijDV/V(it) 

t0.lM,.N/CoNiTj/P I ,T«OhI ,PflOlAu,FHU,uEtRtE ,5AR ,CAR ,□« I ,NOlH 
CPhHoN/HtA) C/AI5,ATM,AID,AICi(j)iATC^(3l 
ulMFivSlDn tVfcOo( J,3) ,F<EiO(3 ,31 , tYtO 1 ( j , 3 ) , EtE U ( 3 , 3 ) , 
ULUlll,3l,cLnn.3> 

i/ItiEi,Sin., P ( 3,i) ,Hl 1 ,3) ,tYtO(3, 3) itYtl (3,3) ,tYt2(3,3» ,EYE3(3,3» 
uI,(EuSiO.J uUl ,3 ),OTl(;>, 1 ) , AOl I 3,1) ,Alo( 1 ,j) ,A 11 n , 1 ) ,AA(A,A) , 
ImIXO 13, 1 ) ,cYi:A( 3 ,31 , 11(3,11 ,12(3,1 ) , 13 ( 3 , 1 1 , TA ( 3 , 1 1 ,15.(3 , 1 1 , 
216(3,1), 17 (3, II, 18(3, l),19l3,l),llu(j,l), 111(3, 1),1U(3,1), 
311^(3, 1 ),TaA (3, 11,113(3,1), 116 (3, IKL.l'ECKi, 3), 

4nHU( 3 ) ,Rn31 ( 3,1 ) ,RiiDi( 1,3) ,Rh0t>(3,3) ,02(3), GiDUl (3,1) 

Ul ME AS) Cm L015I (3,1 ) ,E11S1(3,1 ) ,tl 3 l (1,1 ) ,aA1(h,<,) , 

IlLiAI ,rtlMA) ,tSrAK(<,,i),lLRJUt(<i,l) 
uIMEnSlCfv V ( 11 ,C(:RY ( 1) ,Al3,3) ,D ( 3 , 3 ) ,5(3,31 ,uKt.O (3 1 , w( ol 
ulHEnSlPri aC(3) 

lOuI vAlEWCE (aCi ,aC( 1 )) , ( Af2,AC( 2)) , (XC3 , aC( 31) 
t,lMEi\S10.« (15(3), OH(3), 0013) 

cOuIvALFACb (uS( 1 ) ,LP() , ((.5(2 ),5P1 1 , (U5( 3) ,ZEkD 
tPi,I»ALF.,CE(jMl 1 ), 31^1 ,(t*'( 2), 5 IjC 12), (UM( 3) ,011112) 
cOulvALEuCtdjDll i ,3P2 CuL ), (U0(2) ,CP2CDL1.,(UD13).m 3DLK) 
ulMEw5in,l OS(3),».H(3),(.D(3) 

LiInEh51D>« Y5(3),YM(3),YD(3) 

i. 

Al = ( 1.0-tP5**2l 
a2 = Ai‘>»1 .5 
a3=Ai*»3 

aA= ( 1 .0+tP5eC05 ( Y( 121)) 

a 5 = A ‘,**2 

a6 = Ah*«'3 

A7=PrtEo A*A5 /a2 

(,= (nHEl»A**2 l*At>/A3 

K5 = Ai*Al /AA 

uUrt = -5uR 1( 1 .C-( EKRUS/R3 ) «*2 ) 

(. 

(.AlL FvlLtR(Y,A,B,5) 
lSoN=l 

5lvE3=3AK*U(3,l) 

1F(SIGE3.G1 .(,U/,) 15UN = 2 
4.MtG(l )=Y(1 ) 

' lMcCiZ) =Y( 2) 
b‘'t:G(3)=Yl3) 

GFl0u 1=Y(A 1 
uFi-=Y(5) 

(. 

PH)1=^V( 1u)»RaD1Ai\ 

PHi2=V( 11) *RADiAH 
lMklAlsV(l2)»RADlAn 
lHtlA2 = V(l3)*RADlA!Y 
i/lj 1 K = l,3 

' 1 aC(K) = 3Ak*A(i,K)+CaR<-A(3,K) 

i.ll = (.P5(lHtlAl) 

(.12=005 ((lHklA2) 

5ll = 5IW( (HtlAl) 

5l2=5li»l (HtlA2) 

LP)=L03 (PHll.) 

LP2=C03(PH)2) 


iPl=iliM(PHll > 
iPi=Sl«(PHJ2) 
CDl.=tOS(OEu) 
iDL = i I N ( uEl ) 
iTlC12 = -iTl»(.T<l 
l-IUCi2 = CU*CT2 
SPiCbL = SP2*Ci,L 
tP2CuL=-CP2*L0L 
iDL6i=-iOL 
iEKf = 0 .0 

iPdiUN.cO.l )Cu iO 7 


X 

lAUlLA-Ic 

^UH 


ON 

4jOjMS 

FOK 

iHtRrtAL JEFOKMATiGN 

c 

XI5 

AS 

TriF 

SuN 

angle 

ON 

Sli)E 

BOOM 

IN 

RADiANS 

c 

aIM 

IS 

THE 

SUN 

ANGLE 

ON 

MAIN 

BdOH 

IN 

RAniA.^5 

c 

AlU 

AS 

The 

SUN 

ANGLE 

On 

LAHPtR BUOM 

IN RADIANS 


CPiXlS=CPl»XCl+SHl*XC2 

Xli=ARC0i(t,CiSX15) 

-01iXji1=SI2*Xa-Sll'>CT2'»XC2 + CTUCT2*X(,3 
aIM^ARCO^KC^XIH) 

CnsxiO = SP2»CiJL*XU-CP2*CLL*XC2-SUL*XC3 
AlO=ARCQ^(CO^XiDl 

XTi = ATCl (U*SIN(Xls)9( i .-ATC2I1 )*CUSXli ) 

YTm = aTC1(2>*SIiJ( AlH) «( 1 .-ATC2(c )*CCj$X1M ) 

tlii=ATLH2)*ili^iA.lUl»(-i.--ATCf;( »CuSXILi ) 

CALL VfcCTRlUC.Ui.GSI 

CALL VECTRl (XC»Un,GH) 

call VECTRL (XC.UL.bD) 

iFlNblM.EO. li)ur TU 5 

uFKYnj) = (tTS-Y(i3) l*0«l 

lEjKYII^] =(tTil-tlWJJ90h.l 

uERYlli)=(YTij-Y(13) )*ORI 

uD 4 K=l,3 

YSlK)=Y(13)«uSlK) 

tM(Kl*Y(14)*oMiKl 

YD(K)=y(15)*0D(KJ 



CO Tu « 

5 LiO 6 K = l,3 

YSlKL=YIi*aS (K1 

yM(K)=yTrt*0MlK) 

yD(K)=YTj*oD(KJ 

b CPhTlMUF 

00 Tu 8 
C 

U>LJ.,xI=VT4L«lP.L- . 

HU,<;)=V(4)*SPl 

P(1,:>)=0.0 

P(t.i)=V(5 )»AT^ 

P(<:,^)=-V(5)*ST1'»CT2 

P(2,8)=V(5»*CT1*CT^ 

-Oll.^ljLi V 

8 t'(i,i)=v(4)*cPi*ys(i) 

P(it8)=Vl4)*iPi+ySl2) 

Kl.:i)=Yb(i) 

c 

P(2,l)=V(5)*5T^; + yM(l) 

£12 .2J_= - V15J *5 1L2C12+.YJ1 i c) . 

P(2,J)=V(5 l*CTi*CT2+YH(3) 
y P(3,i)=0.0 
t'(3.^)=0.0 
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H(3, j)=-V<6» 

C 

iJU.L_ LttEitToi H X 

L 

LMO = V(l )«V(2) + \<(3) 

cJil=XU L 

tMfcAK = £MO'»tH l/( EMO + EHl ) 

L 

oV, liO 

no cLol n , n=-H (1 , 1 ) 

1F( liUw.tO.DC.U 10 112 

ELil il ,1 l=y l«.ljt»SiL2»CLLL-Y.i:Li.l . 

LLil 11 ,2)=-V(A)*CPc*COL-^D(2 J 
£U1 (1 ,3J=-V(4>*5DL-yD( 3J 

uE Tu llj _ - - 

il^ cLil(l,l)=V('.)«SH2«Cl.l. 
tUl (1 ,21=-V(A)stPi»CDL 

Hi 11 n .3 1 = -mAL»MlL 

C 

llJ call I«EKT(EM8AR,ELCil,cLCl,EYE00» 

LAlL- I«*£-KT(--cMoAK,tLll,ELOl ,EVc01) 

lALL IhEKTl-tHbAK.tLCl.ELll.cYElU) 
lAlL lilERTlErtPAR.EUll.fcLil ,EYEin 

I/O llS 1=1,3 
uD li5 J=l,3 

c.YEJ.lL^JJ-=tYEOtl,J»+cYtOLll, J) + EYE01(1. J»+tYtlOU,J)+EYElltl,J» 

tYE2lI , Ji = kYtOi ( 1, J)+EYE1U I , J) 
cYe3(l,J) = tYtU(i,J)+ElEiOlI,J) 

lib cYtAll .Jl = bYcQil.t.JJ*LYlLaCLl.altEYEllUtJl _ 

C 

ol( 1 ,11 =LP,; 

kJ.ll,21=3PA 

old ,3)=0.U 
C 

iJ_La(LI_=1i.1 

LMcGi(i,n=OMEi,(n 
12U bTl(l,i)=Gi(l,l) 

_i 

c 

L i-AlClLATE 01 uCl 13 = Alo aT 122 

L U3 3Xi_ __1a3_ 

C 

uC 1^2 1=1,3 

AlUti,U_=0.0 

oO lc2 J=l,3 

I2i Alu(i,n=A10ll,n + ol ll,J)«EYb3U,ll 
C lAlCvLaTl 12 DuT GTl = API A l_ jT AT EP t NJ„1 2 1 _ 

C 3a3 jXi 3X1 

t calculate 111 uOI uTi = AilJ *T STATEMENT d3 

uT. .1.4.3 1. = 1 ,_3 

hOi(i,i)=6.6 
Alio 1 I , 1 )=U .L 

tiii-1^3 j=u.a, .... . . 

l2i A0Hl,l)=A01lI,l)*cYc2lI,J)*cTl{ J,l) 
i23 MllO(I,l)=HliOl I.ldcYElildJloGTlU.l) 

_t .C.A.LC.U.LATC _.0.I j 2 Cr Alio = All AT 12** 

Mil (1,1) =0.0 
uH 1,:4 J = l,3 

i2A Alll..ll=Ailll .ldulll.Jl«'AliClJ.ll . 
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L 


uO 1^5 1=1,3 



AA( I ,J) =tYEl (1 , J) 
AAtIfA)=AOi(l,i) 


12& AA(4,A)=All(l,n 
C 

mMU4-IX^1^34 * 

KHd(2)=S(2,3) 
kHu( 3)=S13,3) ^ 

L 

uO 130 1=1,3 
KHU1(1,I )=AHUm 

U=^3^>ilU^aU4- - 

X3u 02(1 )=01(1,I ) 

L 

UU.-L T 1 

CALL VtCTHDrtE0,£Yt4,T4i) 

CALL VtCTllG^,tYEli,T31 

call Vfc-Cl^lELai,UHtG^OWEO-,-cLlUTA) 

CAlL VECT2(EL0i ,*-HtG,G^,tLll ,T5) 
call VECr2(EL0i,o2,ChEb,tLll,T6) 

XJU>L-,ViX.Ui£LQI.»A^2-».(U , £i-4l-, 1 11- 

call VECT2(ELll,LMEG,0rtEC,EL01,Ttt) 

C 

L- . - cAlCuLATL cvHuT OUT RhCi = KHUD AT STATEMENT 132 

V. jXi 1 a3 jXo 

uO 1 j 2 i = i , ^ 

uILlJi J=i^ — - - 

132 kHuD(I , J)=KHuT( I ,n«KHul (1 ,J) 

C 

013 U5 J=1,3- . - _ _ 

x33 aH^- 0 IJ, J ) = i . U + kHlD ( J , J ) 

c 

iALl-_V(LC_I3lEiJ}li<vHu(Li.ELUiJ9J _ __ _ 

(.ALL VEC I3(ELli ,KHi.0,EL0i,Tl'0) 

LALL VEC rBlDhEol ,EYEil ,Gi, ri j) 

LAlL VtCT31Gi,EVcli,LHcCl.Tl‘*) 

C 

LAlL VcC II (RnO.EYEll ,T15 ) 

lAlL VE CTUCKEb.cYcn .Il t) 

L 

wlOPKl ,l)=-5Pi*LMEG(3) 

GlurKii.Xlr.CKZ^rMEGl 

LluOrt j,l)=SP2»rhEo m-CP2* uHtGI2) 

L 

LA i..Cv.LA.Ti: Lij D.UlJU.DU_.5..JLiJ ftSSU- 

Ill OlT GIDlT = T12 BtLuW TJ ilATtHENT 136 

uP 1 j6 1=1,3 

lli(l,l)=0.0 

1 K ( I , 1 1 =Q.C 
i.0 136 J = l,3 

.HL(.iiiL=Til(IjU>cYEOi{J_,Jl‘>G‘lDL.TUin _ 

111 ( 1, 1 ) =Ti2 II , 1 ) + cYtli ( 1, J)?G10UTU ,1) 
i36 lOmT iNuE 

V 

./P IaO 1=1,3 

tOiS n i , i) = 3 .0*G*Ti ( 1 , i)-T 2 ( 1,1 )-DELJn* 1T13(I , n + Tl<*( 1 , i) ) 

J -J L.E UjXI tEfl par ■>T^ U ii L* Erf&KR AQtUiQJAlIS Qi 1 J _ 
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2 + T6I i ,i) ) + t^'fcA^*l06LwDT#«2)■> (7l 1 ,1 McHBftk* 18 II »1 )+tK6AR*o 
3<-(T9l i ,i) )-uE.-OlT*( TIK 1 ,1 I+T12U f in 

iAU tliSJ U ,i)=3 ,0*t»TA5 11 ,n-Tlbll ,1 ) -DfcLOO I* I T i3 U < il + T 1*. U . 1 IJ 

1- (UElD[jT»*i' )*TjU ,1 ) ♦EnBAR»To( 1 , U+EheAR»&*TiOl I ,1) 

2- DcLum*n2( I.l ) 

C CALCcUTt 01 PtT FiliT * ElSl AT STATEMENT lAA 

«. U3 3A1 1X1 

-cliTU, 1 )=u.o 

wP 1*,A' J»1 ,3 

lA*, tliTU,n = tloTll,n+v.l(l,J)ecllST(J,l) 

<4 >14-1 ♦44--fe-UJ 1-W1-) »-V-< 1-84 syf-LC*.I-=^V 1 1 !>4-*-0tL 

«. 

l/P 1 h5 1=1,4 

4C -ms J = 1 ,4 . . 

i4> mAI t 1 , Jl =AA( 1 ,J) 

V. 

uH la0_.lja^3 

i5y tSTAN(l,i)=ECliTll ,1) 
tSlAKl<,,l)=ElSr(l,l) 

LAlL IrtVtRT ( aAI ,4) 

(. 

1 LA.CLitTf AA I Ou T F:,TAR _s_lii£LU IJfc AT ^ T AT fMrNT iS.J 

C 4X0 hXI 4X1 

uP l5l 1=1,4 

JDkflxitlJ ♦41=0.4 - 

on li.1 J = l,4 

l5l lOKOUE (1,1 ) = 10 k0LE ( I , 1 ) ♦hA 1 ( l , J ) »E AT AR 1 J , 1 ) 


'«(l)=Atl,2»*uKtGl'3)-Atl,;)*tlME0(-i) | 
)=-A(l,l)*nrtEG('4)'*A(l ,3>*bMcGm 

-liJail^ilo.LKiCaai-AiX.i) *DM£b-Ul 

«(4 >=A(2,2)*L'KcG(3)-A(^,j)*0rtEo(,!) 
b( 5) =-a( cl ) »DrtEG( 3)*A( 2 ,3) *uMtG(l» 

ri(6 )= A(2 .1 )»bW^G(21-A(^. ,;)*P HEoll) 

C 

bO 2G0 J=l,4 

♦:0u 0 EnY ( Jl = tOKOLE( J.l ) 

oEnVISi =V(0) 

V. 

oP 2iD J = 1.6 

^lo uEk Y 1 J + 5 ) =u ( J ) 
oERVili )=A7 

’RcTORu 

EIlD 
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r» r% r» 


iUdRLUT I^B IhPuT (DEPfcND ) 
i 

^.0^^'^.^/BLDLl<l/A^,B^,AZ,C^El.A,EPS,N^l,lSU^l,ALP^A,ERRCS 
CCPULN/CLNiTi/k*! ,T»CEI ,P#DlAA,Ff<l,,uE&REE ,SAR ,CaR 
'v.OMMLN/X»N^/ALEAt fdETAt ,LAMAEfPSI,LKDCl3) 

Jlt'EASlOh t)EPEhO(9>,iXLMC(3), PB(3,3) 

OlMEASlON i ( 3, J) iSfl( 3 ,3) , CPC13),CPEC(3) 

C 

ALEAt'ALFAE»RAUlAN 
bETAt=BETAEeRAuIAN 
cAHAt-GAPAE«RAb IAN 
C 

iALFAE=SlN( ALF aE ) 

CAlFAE=COS I ALFaE ) 
iBETAE = SlNLBkTAE) 

C8tTAE=CuS(EETAEl 
iGAPAE^SiN IGAMAE 1 
GCAMAE=CCStCAHAE ) . 

C 

ANGl =CBETAE*GGAPAE 
aN.g 2 = CBEUE*SGAPAE 
ANC3=SBETAEe3GAPAE 
ANGA=SoETAt*CGAPAE 
C _ 

SU,i)=AhGl + iALFAE«AhG3 
i(2f 1 )=-ANG2+SALFA£*ANGA 
i(3,i)=CALFAEe3BETAE ■ 
b( 3 td )=CALFAE«:>CAHAE. 
i<2W)=CALFAt*CGAHAE 
3(3i2)=-iALFAE 
5( 1 13 )=-ANGA-*SALFAE»ANG2 
3(b,3)=AAG3*3ALFAE»AAGi 
. . , 5(3,3)=CALFAE#iBETAE 

C ' 

PSl = PSI*RAi)IAN 
LPiI=CUS (PS_IJ 
iPSi=siN(p'si ) ■ 

C 

6B13,n=CP3I 

BB13,2)=SPSI 

oB(3.3t=C.O 

tsBlJ ,n=-BB(3t2) 

6B( i ,2)=BB('3 ,1 >' 

B6(1,3U0.0 

..6B(2,1 )=0.0_ . 

bB(2,2t=C.O 
BB(2.3)=1.C 

rAkCLLAlE ~BBT obr'sT = 5A AT STATEMENT 31 
3X3 3X3 3X3 

DC 31 l = ... 

DO 31 J=l,3 
SAU.JMC.O 

. DO. 3.1 K=.l,3 

31 SA(1,J)=SA(1 ,JMBBrK,I)»5(J,Kr 
C 

Al = ( l.rEPS.««2> 

A2=(1.-*EPS*CPS1 ) , 

A3=SwRT lAl » 

aa=cfsi+eps 
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k'tD^^.S~Ai>^Al/A^ 

kRX:RAi>IuS«CPSl 

KRt=KADlLS*SPSl 

VX=-«Z*OPEOA*5P51/A3 

VY = A4*IjH£&A*AA/A3 

lHC«1)=0.0 

ut»Cli) = (RRX*»y-RKY«VX)yRA01Ui®*2 
LRCCj)=0.0 

1 ;D 50 lclf 3 
3XuHC(l )=0.0 
ur 5C J=l,3 

50 ^XuHUl )=5XDHC( 1 >*>t J, J)«0HC ( J> 
fO 70 1=1, J 

70 uHE01I)=l,HBCm'»3XLfCI 1)«DEGREE 
i~U 

L »0 8 L I = 1 » ^ 

I'C 80 0 = lf3 

L=L + 1 

8U uEPE^D(U = SA(I • J) 
uC 9C i=l»i 

OHfcG(n=LMEC(n«RAuI/iN 

90 uEPEKDl I+6) = LMtCH ) 

HETUKN 

tNi) 


r' o o o 



iMER TO* «••*•«•***•* 

C***««o***'«i»*««»««olNE»TlA DTAtIC IOTA ZERO********** 

C 

5UaRLUTlA£-.lA£RTC(P.htET£0) • • 

CCKHCN/BCOX/'»Ut> . ' 

UICEASION P(3,3),HU,3»,EYE0(3,3l,B(l,3» ,S(1,3», ' ' 

»RU3jJJjSm,3J,lilJ.3J _ I •.<: -■ 


CAA.CUATE tr. ,l!0T P = h AT STATEHEAT 5 , 

1X3 3X3 1X3 

DO 5 I=lf3 
h(l,l)=0.0 
00 5 J=lf3 

5 Hl,i)=HU, I ) + V( j)*P(J, 1 ) 
tHM = vU)*Vl2)tYr3) 

EHHl=I./EHh 
uG 1C l=lt3 
H(l,ll?»lll-.1J?EMPI 
10 COATiNUE 

c 

00 20 -J = i.3 ,1 

B(1,J)=P(1. J)-H(1,J) 

3(i,j)=P(2>J)-rl(i»J) 

T(1,jI=P(3,J)-H(1,J) 

C 

20 CONTINUE 
C 

C CAlCcLATE iRT CCT B)«EMlti) =«1 AAO 

C JST DuT i)*EMl,2) =5 AND 

C ITT OCT l)#EP(lf3) =T1 

C 3X1 1X3 jX3 below TC STATEHEAT 3C 

CO 30 l=lf3 
DO 3C J=lf3 
C 

Kl|I,JI= Rll,I)»R(l.J)«V(l) 

311IiJ)‘ S( 1 ,nes(l ,.)»V(2) 

T1(I,J»* T(l,I)*l(l,j)*V(31 

0 

tYEOd >U)=K1 (I .J)*X1(1 ,J)«T1 (I, J) 
tY£OiI.J)=-ETEo(l,J) 

3u COwTiNuE 
t 

c=(R(l ,1)**2-*R( 1 ,2»**2-tR(lt3»*'»2J*W(l)* 

1 (SI 1,1 »**2 + Sn, 2) *#2«S( 1,3»**2 »•«(«:) ♦ 

2(T( 1 ,1 J**2*T (1 ,2 )**2'»T( 1 ,3 )**2)*V( j ) 

C 

JO AC 0=1,3 

tYEC(J,J) = EYEO(J,J»*do*t»+U 
AO CONTINUE 
wETUkN 
tNU 
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o^5u8i<0UT |N£ 

Ctooo0o*«ofto©o«o« o^CuFl^UTEi ELEMENTS CF INERTIA JYAO IC 
SUbRLUTiNE iNEKTtEMBAR.ElA.ELB.EYE) 
tlPENSlCN cLfl<l»3),ELB(l,3),EYt(>»3),U(n 
C 
c 

C LALCLLATE IElA TkANSFOSE OoT ELBI'^Ef^bAR = EYE AT 1C 
jD 10 l=l,J 

uO lu J-I 9 J 

-lu tYEU,JJ=-£^^fcAii«(ELA(l , li^ELBCltJ)) 

C 

C CAlCLLATc ELA JCI tLd TRANSPOSE = U AT 1 5 

C 

Ul )*0.0 
dC 15 J=i,3 

15 uU)=Uil)^bLAU»j)«Et8U»J) 

Ui )=U( 1 )«EPdAA 
uC 2C 4=1,3 

2u tYE(w,u)=EtElJ,J)+u(l) 

KETUHN 

ENO 
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(.eeoe«ee»«c««««o««cSu6i>CLT INc VECTl**®** 

t****#» VEtTwR PRCtbcf A»e*fi»7 A ANb I ARE JECTtRi.o iS A CItACIC 
C***9*« X IS VtCTCR PRLOcCIt * 15 bOl PRuDoCI 
SUbRcUTlAE VtCTKArPiZ) 

blHEASlCA A(3),e(3,3),^(2,l) , AAX ( 3 , 3 i , i/UP ( J 1 1 ) 

kC 1C t=1.3 
10 AAx(l,n=:0.C 
AAXI lt4:)=-M(3) 

AAxl 1,9)=A( 2) 
mAX (^,1 )=A (3 ) 

AAx(<;,3)=-Am 
' «Ax( 3 ,n=-A(i) 
mA X ( 3 f I =A ( 1 ) 

C 

c 

C cAlCcLATc B Dbl A = DIM AJ»D 

C 3X3 3Xi 3>1 

L AAX CuT DoK = 2 BELbW AT STATcPEMS 1 a AAC. 15 

I.C lA 1=1, j 
bUMI 1 ,i)=0.0 
bC lA J=if3 

lA bUI1(l,i)=DbP(I,l)'» B(I>J)* AU) 
uO IS 1=1, 3 

X (1,1) =0.0 

OC IS J=l,3 

IS 2(I,1)=Z(I,1)«aAX(1,.;)*CLH(J,1) 

c 

RETURN 

END 
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0 «su8nDL T I Nt ^ECT2 ao*o^« 

C***«#0 xECTOR PkOU’UCT AXBXCXC=Z A,B,C,0 and i ABE VECTORS 
C****** X IS VECICr PREDUCI 

SUeRtUTUF V ECT 2 I A , B ,C , C ,Z ) 

LlHEA'SiCA An,j),BI3)>CI.),D(l,3),AX(3.3)>BX(3,3),t.X(3,3), 
1 DUHll3.1)«OUR2(3.1),Z(3fl) 

t 

uO 10 1=1,3 
AX( I ,I )=0.0 
bX( I ,1) =C.U 
10 m.XU,I)=C.0 

AX( 1 ,2> = -AU ,31 
AX(1.3)=A(1,,:) 

AX(2,1 l=A(l ,3) 

AX12,3»=-A(1,1> 

AX13,1)=-A11,2) 

AX t 3,2)=A( 1 , 1 ) 
t 

dX11,21=-B131 
' dXU ,31=e(2 ) 

,i>xi2,i )=;Dr3) 
oX(2',3)=l-B( 1 1 
bX(3,,n=-P( 2.1 
DX13’,2)=Bm 

c 

».Xl 1 ,2)=-Cl31 
LXn,3)=CI2) 

CX12,l)=C(j) 

(.X(2,3 )=-C (1 I 
«.X(3,1)=-C(2) 
tXi3,2)=Ui) 

C 

t 


V. 

CAuCUaTE 

cx 

Dul 

UT = DUHi 

A 1^0 

C 


Ba 

OuT 

0UM= COH^ 

AnO 



AX 

DLT 

0L(*»= Z 

BtLCJW 





3X1 3X1 



uD 2C 1=1, A 
t,UKl(l,l) = O.C 
bO 2C J=i,j 

20 bUni n,ii=L.U'iii ,1 Hcxd ,j)»u(i,j) 

L.n 3 l 1 = 1,3 
JUH2ll,l)=0.C 
lD ?C J=1,3 

30 i.L>h2(1 , 1 l=uUl<2< I ,1 dEX( I .J)»CUM1< J, 1 1 

liC Au t =1 , j 
i.( 1 , i) =0.0 

00 Au J=i,J 

Ao ll l,i)=Zll,l >4riXtI,J)»LLF2( J,l) 
kETUkR 
EOfll 


l.o*6os#6##so*«s«»««-iteRCl;nhE VECT3 o*os*«** 

Coeoes* VECTOR PRCCUCT »XB«C=Z A,C AND I ARE VECTCRi.B IS A OtAOIC 
C*o*oe* X IS VtCTCk PRODUCT, * li UPT PRODUCT 
SUBRlUTINE VECT3(A,B,C,Z) 

LIHEaSIOA a { 1 , J ) ,B( 3 ,3 1 ,C( 1 , 3 1 , AX( 3 , J) , OCR t 3 il > , X( 3 ,’l ) 


OO 10 1=1,3 
10 AX( I ,I )=C.O 

AXd ,21=-All ,3) 
AX(l,3)=A(l,i) 
AX(2,1 )=A(1 
AXl2,3)=-All,n 
AXl 3,1)=-A( 1 ,2) 
aX(3,2)=A(1,1) 

C 


C 


e 

CAlCLLATE b 

OlT 

CC 

= uUM 

L 


CCT 

CLM 

= L 

C 

3X:i 


3X1 

3X1 


UC 20 1*1,3 
OUHlil ,1 1*0.0 
00 2t J*l,3 

20 uUMlU,i)=DURi(I,l) + b(l,j)« C(1 
bO 30 1=1,3 
X( 1,11=0.0 


AND (CC=C TRANSPOSE) 
BELGk 


J 1 


00 3C J*l,3 

30 i(l,l)=Z(I,l)'»AXlI,Jl»DOMJ,l) 
kETUHN 
tNU 
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^^600* ocSuBKDtT INk 

t^«><=o«<.oc5c«?co soseoM^TKl A INVF^SICK oV GAUiS- JCRl/^^ cL iM 1 NAT I CN*®*’!'® » 
^)UbPlUT I »\E I NVtRT ( A ,A ) 
blHE^5iC^ ,B(5 l) .C (50) tLi.(50) 

jUH s i . 0 
lP 5 1 = 1 ,K' 

5 :»l;MriUhoA( 1 • 1 ) 

C 

kavg = i^..o*«(-alCoiu( sum ) / N ) 

1,0 6 1 = 1, N 
lO 6 J = 1 ,N 

b «( 1 , j)=All , Jl^nAVG 

V 

Ur IL J= 1»N 
1 J LZ ( J ) =J 

0 C 2 u I = 1 , w 
^=1 

V = A ( 1 , 1 ) 

LM-1 
lP = I -t 1 

IFIN-Lk) lA.il.ll 

!i or ij j=lP,n 
» =A ( 1 , J ) 

1F( AtS(Wl-AF^m ) 13,13,iZ 
1 <■ IV = J 

» = X 

1 J lPi\T 1 NuE 

l4 iF(V.LT.l.k-a) Gl to 2oC 
IM = 1 .C/Y 
lC 13 J=ll^ 

G( Jl=AlJ,K) 
w(,j,X)=AtJ, I ) 

*<( J , 1 ) =-!.( J ) »Yl 
Alt , J) =A(1 , J )*YI 
13 c( J)=Al I,J1 
«(1 , 1 ) -Y1 

j = LZ I n 

lZ t I ) = LZ (K 1 
lZ IK l = J 
lO !■, K = 1,N 
iFII-K) i6,1S,i 6 
lb bC lo J=l,n 

IFlI-J) 17,1 cl 7 
17 «(K,j)=AlK,J)-btj)eC(K) 
lo oCATiNoE 
ly LOkTINuE 
2b bCMlKoE 

uP 2b0 1=1, F 

iFll-LZtni lOO.oOO.lOo 
iOu K=l+, 

bP 5l0 J=K,K 
iFll-U(j)) tOb,oOb,iCO 
cOO F=l.ZlI 1 

lZ I n = LZ U) 

lZ tJ ) =M 
bP 7v,0 L=1 ,F 
i.(Ll=At I ,L) 

A( l , lI = A iJ,L ) 
iOj h(^,l)=CIL 1 
5^0 V-PuTlFluE 



iOO ^.C^TINuE 

L MAKfc IT A SYMMETRIC MATRIX 

C;0 2S0 I *1 ,K 
l/C ZiO J = I,N 

AVO=lA(I ,J» + AU,n l*.5*RAVl, 
All ,.)=AVC 
h( J , 1) =AVC 
C50 CONTINUE 
kETURN 

c 

RbO N = -1ABS(M 
KEIUkN 
tNO 


:>UbRLUTUE ELLkR (Y.A.B.S ) 

C 

1.1 HENS 10l« Y ( 1) ,A (3,3) ,b( 3,3) ,S( 3.31 
t 

B(3,i)=CLS(Y(li)) 

b(3,,;)=5IN(Y112)) 

B(3,3)=0.0 
B(l,l)=-B(i,i) 
b(l,2)=6(3, 1 ) 

B(l,3)=C.O 

P(2,l)=C.O 

B(2,2)=0.0 

cti,3)=1.0 

AU,i)=Y(6) 

A(i,i)=Y(7) 

A( 1,3)=Y(8) 

A(i,i)=Y(9) 

Ali,i)=Yno) 

A(,;,j)=Y(n) 

A(3,i)=Yl7)*Y(ll )-Yie)*Y(lu) 

A( j,*)=Y(8)«Y(9)-Y(6)*Y( 11) 
A(J,j)=Y(6)*t(iC)-Y(7)»Y(9) 

1. 

C 

t. CALCUATt AT DCT bT = S AT 180 

C JX3 3X3 3X3 

L 

tO IdO 1=1,3 
iJC loO J = l,3 
S( 1 , J) = 0 .0 
UC IbO K=l,3 

180 S(l , j)=S(l , J)+ AlK,nABt.,K) 

0 

' KE IUKN 
tNO 



6UBRCUT114E VECTRK AtEtC) 

UlPENSlCh All), Bin, OU) 

C 

C i/ECTRI calculates THE \iECTUR TRIPLE CR OS 5-PRCDUC T- ( AXB ) XB AMU 

L PUTS ToE RESULT IN D WHERE (AXB),XB = (-BDCTB)A ♦ (BOCTA)B, ThE 

C result, D, is NCRHALIZEC— THUS D IS ORThONCfcMAL 

C 

B12 = B(1 )«B(1 ) 
o22=b(2)«B(2) 
b32=fa(3 )*B(3 ) 

Alt)l=Am*B(l) , 

A2B2=A(2)*B(i) 

A3B3=AI3)*6I3) 

l( 1 )--a( 1) «( B22'»B32)-»B( 1 )»( A2B2'»A3B3 ) 

L(^ )=-AU)*lBU+B32 )-*Bl2»*(AlBl+A3b3l 
t(j)=-A(3)*(fcl2 + t22)-»B(3)*(AlBl*A2B2l 
OMl=i.O/(Sum(‘J( l)«D(l)*i:(i)»D(2)*U(3)^D(3) )) 

D(1)=D|1 )*OPI 

U(2»=0I2)*0P1 

D( j)=D(3»*DMl 

RETURN 

END 
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Program Output 



For each set of input data the run number and all of the initial 
primary satellite parameters are first printed out.' The variable parameters 
are then computed and subsequently printed out at 2- minute time intervals 
throughout the 440 minutes of simulated satellite orbital motion. The com- 
puted printout includes: 


TIME, min 
R, km 

RX, km 

RY, km 

ROLL, deg 

PITCH, deg 

YAW, deg 

ECLIPSED, 
dimensionles s 


Time in orbit after specification 
of initial values. 

Magnitude of the radius vector 
from the earth' s center to the 
satellite' s center of mass. 

Rcos ijj 

Rsinijj 

Roll angle. 

Pitch angle. 

Yaw angle. 

Occultation parameter. 
ECLIPSED = 2; satellite is in 
sunlight. 

ECLIPSED = 1; satellite is in 
earth' s shadow. 


A typical printout follows. Only a small portion of the computed 
time history is shown since the remaining output has the same format. 


RON 1«»9 


Ml, SIDE BOJ'* AND MA,UAMPER BOOM TIP 

HASS 

= 

0.1175A 

SLUGS 

m2,maih booh tip mass 



0.1592 

SLUGS 

M3, PAYLOAD MASS 


= 

8.7999 

SLUGS 

LI, SIDE BQuM AND LA, DAMPER BOOM LENGTH 

= 

35.00 

FEET 

L2,MAIN BOOH LENGTH 


= 

60.00 

FEET 

L3, DISTANCE FROM 3AYL3AC C.G. TO HINGE 


1.375 

FEET 

01, SIDE BOuM AND DA, DAMPER BOOM DIAMETER 

= 

0.25 

INCHES 

D2,maih boum diameter 


= 

0.5C 

INCHES 

HI, SIDE, AND HA, DAMPER, BOOM THICKNESS 

= 

1.40E-03 

INCHES 

H4,HAIN B3uM TilCK'^ESS 


= 

2.00E-03 

INCHES 

Ui.l) •PAYLOAD ROLL INERTIA 


=. 

4.00 

SLUG-FT2 

1(^,2) fPAYLOAO PITCH INERTIA 


= 

4.0C 

5LUG-FT2 

fPAYLOAO YAW INERTIA 


= 

4.00 

SLUG-FT2 

PHilfYAW ROTATION 3F SIDE BOOM 


=. 

-30.00 

OEG 

PH12,YAW ROTATION OF MINOE AXIS 


= 

120.00 

DEG 

THETAltRQLL ROTATION OF MAIN BOOM 


= 

0.0 

DEG 

THETAZfPITCH ROTATION OF MAIN BOOM 


= 

0.0 

OEG 

DELTA 2ER0,NULL POSITION OF DAMPER 

BOOM 

= 

0.0 

DEG 

KDtDAHPER SPRING CONSTANT 


= 

85.00E-05 

FT-LB/RAO 

CDtDAMPER DAMPING CONSTANT 


= 

0.395 

FT-LB-5EC/RAD 


ECCENTRICITY 
SErtl-MAJOR AXIS 
EARTH RADIUS 
INITIAL TRUE ANDNA^Y 
GRAVITATIONAL CONSTANT 
OMEGA, MEAN ORBITAL RATE 
SUN ANGLE 


= 0.0 

=7302. A3E 00 KM 
=S37B.lbE 00 KM 
= 0.0 OEG 

=3986. 13E 02 KM3/SEC2 
= Un. 75E-06 RAD/SEC 
= 30.00 DEG 


ABSORPTIVITY 

LINEAR THERMAL COEFFICIENT 3F EXPANSION 
THERMAL CONDUCTIVITY 
HEAT RADIATION OF THE SOURCE 
THERMAL LAG THE 


0.70 


=10A0.00E-08 
=A167.00E-03 
=3065.00E-03 
= 0.0 


IN/ IN-F 
BTU/HR.-IN-F 
BTU/HR-IN2 
MIN 


INITIAL ROLL ANGLE 
INITIAL PITCH ANGLE 
INITIAL YAM ANGLE 


= -10.300330 OEG 
= -30.000300 DEG 
= 30.000300 DEG 
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TIME R RX RY R3LL PITCH YAW ECLIPSED 

0.0 730^.AR 7R02.A3 0.0 -lO.OOE 00 -30.00E 00 30.00E 00 2 

2.00 7302. A3 72Ad.67 BSA.Al -98.80E-01 -29.6AE 00 29.82E 00 2 

A. 00 7302. A3 7088.20 1755.80 -95.00E-01 -28.59E 00 29.30E 00 2 

b.OO 7302. A3 6823.38 2601.35 -88.A3E-01 >26.88E 00 28.51E 00 2 

8.00 7302. A3 6A58.09 3A08.59 -79.09E-01 -2A.58E 00 27.A9E 00 2 

10.00 7302. A3 5997.73 A165.S5 -67.23E-01 -21.75E 00 26.35E 00 2 

12.00 7302. A3 5AA9.07 A861.39 -53.35E-01 -18.A9E 00 25.17E 00 2 

lA.OO 7302. A3 A820.18 5A85.55 -38.16E-01 -1A.88E 00 2A.03E 00 2 

16.00 7302. A3 A120.3A 6028.95 -22.55E-01 -ll.OAE 00 22.99E 00 2 

18.00 7302. A3 3359.83 6A83.59 -75.31E-02 -70.72E-01 22.09t 00 2 

20.00 7302. A3 25A9.86 68A2.78 59.51E-02 -31.30E-01 21.3AE 00 2 

22.00 7302. A3 1702.35 7101.23 17.06E-01 76.11E-02 20.68E 00 2 

2A.00 7302. A3 829.78 7255.13 25.18E-01 AA.OlE-01 20.06E 00 2 

26.00 7302. A3 -55.01 7302.22 29.96E-01 77.22E-01 19.39E 00 2 

28.00 7302. A3 -938.93 72A1.80 31.27E-01 10.6AE 00 18.58E 00 2 

30.00 7302. A3 -1809. lA 707A.78 29.2AE-01 13.09E 00 17.53E 00 2 

32.00 7302.A3 -2652.65 6803.59 2A.11E-01 15.0AE 00 16.19E 00 2 

3A.00 7302. A3 -3A57.il 6A32.25 16.27E-01 16.A7E 00 1A.51E 00 2 

36.00 7302.A3 -A210.67 5966.21 61.32E-02 17.38E 00 12.A6E 00 2 

36.00 7302. A3 -A902.2A 5A12.35 -58.A7E-02 17.79E 00 10.05E 00 2 

AO. 00 7302. A3 -5521. 6A A778.79 -19.21E-01 17.73E 00 73.11E-01 2 

A2.00 7302. A3 -6059.75 A07A.90 -33.A9E-01 17.25E DO A2.66E-01 2 

AA.OO 7302. A3 -6508.66 3311.02 -A8.23E-01 16.A0E 00 96.13E-02 i. 

A6.00 7302. A3 -6861. 7A 2A98.39 -62.95E-01 15.2AE 00 -25.57E-01 2 

A8.00 7302. A3 -7113.81 16A8.98 -77.18E-01 13.8AE 00 -62.A8E-01 1 

50.00 7302. A3 -7261.16 775.29 -86.63E-01 12.25E 00 -99.15E-01 1 

52.00 7302. A3 -7301 .60 -109.80 -9A.80E-01 10.58E 00 -13.69E 00 1 

5A.OO 7302. A3 -723A.56 -993.28 -95.82E-01 87.30E-01 -17.29E 00 1 

56.00 7302. A3 -7061.01 -1862. lA -9A.A5E-01 69.06E-01 -20.95E 00 

58.00 7302. A3 -6783.51 -2703.59 -91.82E-01 51.16E-01 -2A.5AE 00 

60.00 7302. A3 -6A06.15 -3505.23 -85.86E-01 33.A2E-01 -28.06E 00 

62.00 7302. A3 -593A.A3 -A255.2B -81.0AE-01 18.28E-01 -31.57E 00 

6A.00 7302. A3 -5375. A5 -A9A2.68 -73.77E-01 37.26E-02 -3A.89E 00 

66.00 7302.A3 -A737.27 -5557.32 -66.3AE-01 -82.93E-02 -38.10E 00 

68.00 7302. A3 -A029.36 -6090.13 -57.98E-01 -18.70E-01 -Al.lOE 00 

70.00 7302. A3 -3262.13 -6533.29 -A9.71E-01 -26.55E-01 -A3.92E 30 

72.00 7302.A3 -2AA6.88 -6880.28 -A1.A7E-01 -32.32E-01 -A6.51E 00 

7A.OO 7302. A3 -1595.61 -7125.97 -33.86E-01 -35.63E-01 -A8.89E OC 

76.00 7302. A3 -720. 8A -7266.76 -27.10E-01 -36.7AE-01 -51.J3E JO 

78.00 7302. A3 16A.53 -7300.57 -21.59E-01 -35.60E-01 -52.96E 30 

80.00 7302. A3 10A7.A9 -7226.91 -17.59E-01 -32.A3E-01 -5A.66E 00 

82.00 7302. A3 1915.02 -70A6.85 -15.j5E-01 -27.39E-31 -56.17E 00 

8A.00 7302. A3 275A.37 -6763.05 -1A.99E-01 -20.76E-01 -57.A8E 00 

86.00 7302. A3 3553.16 -6379.70 -16.61E-01 -12.8AE-01 -58.63E 00 

88.00 7302. A3 A299.65 -5902. A1 -20.1AE-01 -39.6AE-02 -59.62E 00 
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90.00 

7302.43 

4982.83 

-5338.24 

-25.45E-01 

54.B6E-02 

-60.49E 

00 

92.00 

7302.43 

5592.67 

-4695.47 

-32.32E-01 

15.11E-01 

-61.23E 

00 

9«f.00 

_.7302.J»3 

6izoa7 

r3983^59 

-40.43E-01 

24.52E-01 

-61.88E 

00 

96.00 

7302.43 

6557.57 

-3213.07 

-49.41E-01 

33.31E-01 

-62.45E 

00 

96.00 

7302.43 

6898.43 

-2395.24 

-58.84E-01 

41.09E-01 

-62.94E 

00 

100.00 

7302.43 

7137.73 

-1542.14 

-68.31E-01 

47.51E-01 

-63.36E 

00 

102.00 

7302.43 

7271.96 

-666.34 

-77.42E-01 

52.23E-01 

-63.71E 

00 

104.00 

7302.43 

7299.14 

219.26 

-85.80E-01 

54.96E-01 

-63.99E 

00 

106.00 

7302.43 

7218.86 

1101.63 

-93.16E-01 

55.46E-01 

-64.20E 

00 

108.00 

7302.43 

7032.30 

1967.79 

-99.27E-01 

53.55E-01 

-64.33E 

00 

110.00 

7302.43 

6742.22 

2804.99 

-10.40E 00 

49.12E-01 

-64.37E 

00 

112.00 

7302.43 

6352.88 

3600.88 

-10.72E 00 

42.14E-01 

-64.33E 

00 

114.00 

7302.43 

5870.01 

4343.77 

-10.89E 00 

32.67E-01 

-64.19E 

00 

116.00 

7302.43 

5300.74 

5022.71 

-10.92E 00 

20.87E-01 

-63.95E 

00 

118.00 

7302.43 

4653.42 

5627.71 

-10.82E 00 

69.71E-02 

-63.61E 

00 

120.00 

7302.43 

3937.60 

6149.86 

-10.61E 00 

-B6.89E-02 

-63.17E 

00 

122.00 

7302.43 

3163.62 

6581.47 

-10.31E 00 

-25.72E-01 

-62.63E 

00 

124.00 

7302.43 

2343.46 

6916.19 

-99.37E-01 

-43.67E-01 

-61.99E 

00 

126.00 

7302.43 

1488.59 

7149.09 

-95.34E-01 

-62.06E-01 

-61.25E 

00 

128.00 

7302.43 

611.81 

7276.75 

-91.26E-01 

-80.41E-01 

-60.42E 

00 

130.00 

7302.43 

-273.97 

7297.29 

-87.46E-01 

-98.24E-01 

-59.48E 

00 

132.00 

7302.43 

-1155.72 

7210.39 

-84.25E-01 

-11.51E 00 

-58.44E 

00 

134.00 

7302.43 

-2020.46 

7017.35 

-81.94E-01 

-13.06E 00 

-57.30E ,00 

136.00 

7302.43 

-2855.45 

6721.00 

-80.78E-01 

-14.45E 00 

-56.05E 

00 

138.00 

7302.43 

-3648.40* 

6325.71 

-81.01E-01 

-15.65E 00 

-54.70E 

00 

140.00 

7302.43 

-4 387.66 

5837.29 

-82.77E-01 

-16.64E 00 

-53.24E 

00 

142.00 

7302.43 

-5062.30 

5262.94 

-86.17E-01 

-17.41E 00 

-51.68E 

oc 

144.0U 

7302.43 

-5662.44 

4611.10 

-91.20E-01 

-17.95E 00 

-50.03E 

00 

146.00 

7302.43 

-6179.20 

3891.39 

-97.79E-01 

-18.27E 00 

-48.28E 

00 

146.00 

7302.43 

-6605.00 

3114.39 

-10.57E 00 

-18.36E 00 

-46.45E 

00 

150.00 

7302.43 

-6933.56 

2291.54 

-11.48E 00 

-18.25E 00 

-44.53E 

00 

152.00 

7302.43 

-7160.05 

1434.96 

-12.46E 00 

-17.93E 00 

-42.54E 

00 

154.00 

7302.^3 

-7281.14 

557.24 

-13.15E 00 

-17.34E 00 

-40.54E 

00 

156.00 

7302".43 

-7295.03' 

-320.67 

-13.79E 00 

-16.67E 00 

-38.37E 

00 

158.00 

7302.43 

-7201.53 

-1209.74 

-13.90E 00 

-15.72E 00 

-36.21E 

00 

160.00 

7302.43 

-7002.01 

-2073.01 

-13.78E 00 

-14.82E 0“0 

-33.77E 

00 

162.00 

7302.43 

-6699.41 

-2905.76 

-13.61E 00 

-13.79E 00 

-31.13E 

00 

164.00 

7302.43 

-6298.16 

-3695.72 

f-13.1l£ 00 

-12.71E 00 

-28.32E 

00 

166.00 

7302.43 

-5804.23 

-4431.29 

-12.74E 00 

-11.62E 00 

-25.21E 

00 

168.00 

7302.43 

-5224.84" 

-5101.62 

-12. HE 00 

-10.40c 00 

-22.00E 

00 

170.00 

7302.43 

-4568.52 

-5696.85 

-11.44E 00 

-91.48E-01 

-18.60E 

00 

172.00 

7302.43 

-3844.96 

-6208.20 

-10.61E 00 

-77.93E-01 

-15.10E 

00 

1*74.00 

7302.43 

-3064.79 

-6628.16 

-97.43E-01 

-63.59E-01 

-11.52E 

00 

176.00 

7302.43 

-2239.50 

-6950.55 

-87.96E-01 

-48.11E-01 

-78.98E-31 

178.00 

7302.43 

-1361.24 

-7170.61 

-78.42E-01 

-31.46E-01 

-42.91E-01 

130.00 

7302.43 

-502.64 

-7285.11 

-69.05E-01 

-13.48E-01 

-73.56E-02 . 

182.00 

7302.43 

383.34 

-7292.36 

-60.46E-01 

58.44E-02 

27.09E-01 

184.00 

7302.43 

1263.69 

-7192.26 

-53.14E-01 

26.46E-01 

59.96E-01 

186.00 

7302.43 

2125.44 

-6986.27 

-47.65E-01 

48.18E-01 

90.72E-01 

168.00 

7302.43 

2955.90 

-6677.43 

-44.43E-01 

70.69E-01 

ii.ddE 

DO 
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APPENDIX E 


EXAMINATION OF THE ANOMALOUS BEHAVIOR 
OF THREE GRAVITY GRADIENT SATELLITES 


By . 

Robert L. Goldman 

Research Institute for Advance Studies 
Martin Marietta" Corporation 


(Because this study is required for the clear understanding of the basic 
publication, it as reproduced in its entirety Trom publication TR-71-07c, 
RIAS, Martin Marietta Corporation, March 1971.) 

' , ^ SUMMARY- 

The anomalous oscillatory behavior of three satellites orbited by the 
Naval Research Laboratory has been examined. The satellites, a part of 
the 160 senes of experiments, were all hinged two-body gravity gradient 
configurations passively stabilized about their three principal axes by 
gravity gradient, and damper torques. Two basic behavior patterns were 
observed. ’Either the satellites were stable with attitude perturbations 
less than ±5^ or their behavior tended towards a low frequency rigid body 
oscillation dominated by large yaw motions and in some cases by yaw 
inversions. A causal relationship between sun angle and the character 
of satellite behavior was observed that appears to indicate that thermal 
distortion is a critical factor in a gravity gradient satellite's dynamic 
behavior. The behavior appeared to be further modified by the existence 
of response frequencies that were higher than anticipated values. 

INTRODUCTION 

Threc-axis, passive, gravity gradient stabilization of spacecraft through 
the use of extendable booms has been demonstrated as a practical means 
for providing an carth-pointing equilibrium orientation (ref. 1). The suc- 
cess of these gravity gradient systems, however, has, for certain satellite 
configurations, been inexplicably associated with a low frequency anomalous 
oscillatory behavior. This unpredictable behavior has usually appeared 
as a sustained large-amplitude, rigid-body oscillation modified, in some 
cases, by one or more attitude inversions. A typical example of this type 
of behavior is illustrated in figure 1. Such a performance was clearly 
seen in tlie flight data collected during an initial series of gravity gradient 
experiments conducted by tlie Naval Research Laboratory (NRL) (ref. 2) 
and more recently in Die data collected from the latest series of gravity 
giadient satellites orbited by NRL. 
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This new NRL flight test data is used in the present* report as a basis for 
a farther examination of a passive gravity gradient satellite's low frequency 
behavior. The investigation has been directed towards the collection, dis- 
play, identification, interpretation and evaluation of data from three of 
these satellites, and is oriented to the objective of attempting to ascertain 
the essential ingredients of the behavior mechanism. 



ORBIT NUMBER 

ECLIPSE — YAW PITCH ROLL 


Figure 1 Typical Yaw Inversion, Satellite 163 

The anomalous low-frequency behavior of a passive gravity graaient satellite 
can be viewed in its most general sense as an unstable interaction phenomenon 
involving coupling between internal dynamic properties of the satellite and 
external environmental energy sources. Depending upon the satellite's 
orbital distance and eccentricity^ environmental energy sources (such as 
tliose due to aerodynamics, solar radiation and magnetic fields) may intro- 
duce destabilizing torques that are large when con?>ared to the satellite's 
stabilizing gravity gradient torques (ref. 3). An adequate understanding 
of, the dynamics of the anomalous behavior -is required before any logical 
attempt can be made to eliminate the problem. The task that irises' in the 
present study, therefore, is one that tries to find out which of the many 
internal and external system characteristics clearly dominates the unstable 
interaction phenomenon. Since the flexibility of a gravity gradient 
satellite's booms and the influence of solar pressure and tjhermal bending 
are generally suspected as being principal .offenders in boom instabilities 
(refs. 4 to 9), they have been given principal consideration. 


^ SPACECRAFT CHARACTERISTICS 

The gravity gradient satellites in the NRL 160 series were launched 
together in the latter part of 1969 and successfully placed in a nearly 
circular 500 nautical mile orbit at an inclination of approximately 70® 
to the earth's equator. Their orbital parameters are summarized in 
Table 1. The satellites essentially moved along the same orbital path 
with a spacing of roughly 100 nautical miles between them. The orbital 
period and precession rate were such that the satellites came' close 
(within 5®) to passing over the same point on the Earth's surface every 
14 orbits. 
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TIP WEIGHT (TYPICAL) 




Figure 2 Satellite Geometry, Payloads 161 and 163 Figure 3 Satellite Geometry, Payload 164 


Because of their susceptibility to rigid-body anomalous oscillations and 
their geometric similarity^ three of these satellites^ payloads 161^ 163 
and 164^ shown in figures 2 and 3^ have been singled out for examination. 
These particular satellites were hinged two-body configurations (ref. 10)^ 
passively stabilized about their three principal axes by gravity gradient 
and damper torques. The satellites were asymmetric with respect to both 
their geometrical shape and mass properties^ but each had a plane of , 
geometrical symmetry with respect to their principal axes. They each used 
three long extendable booms with a passive hinge damper attached to one 
or more of the booms. Tip weights were located at the deployed ends of 
the booms. The long booms and tip weights were needed in order to obtain 
a large enough moment of inertia for gravity torques to be effective^ 
while the damper mechanism was designed to dissipate energy in order to 
inhibit tumbling and limit librational motions. Although active devices^ 
such as momentum wheels and thrusters^ were available on these payloads^ 
their use was not required for stability. The important physical char- 
acteristics of the three satellite configurations are summarized in 
Table 2. The inertial properties listed in Table 2 are consistent with 
the definitions used in the formulation of the general equations of motion 
for a hinged two-body satellite given in reference 10. 

Satellites 161 and 163 


Gravity gradient satellites 161 and 163^ illustrated in figure 2^ have 
basically the same geometry. They differ mainly in their boom construc- 
tion^ a dissimilarity that makes examination of their motion attractive^ 
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Table 1^-Orbital Parameters 


Parameter 

Valve ' 

Eccentricity 

Inclination 

Period 

Pengee altitude 
Apogee altitude 
Orbital precession (eastward) 

0.00203 

70.014® 

103.46 min/orbit 
490 5 naut miles 
506.5 naut miles 
2.121 deg/day 


Table 2— Satelbte Physical Properties 


Property 

Units 


Satelhte 

- ‘ 

161 

163 

164 

Payload weight 

lb 

235 

247 

283 

Mam boom tip weight 

lb 

5.28 

5.28 

5.12 

Lateral tip weight (each) 

lb 

3 76 

3 76 

3.78 

Mam* boom weight , 

gm/ft 

6.866 

4.584 

6 866 

Lateral boom weight (each) 

gm/ft 

2.709 

4 584 

2.709 

Mam boom length 

ft ■ 

60 

60 

60 - 

Lateral boom length 

ft 

37 

37 

35 

Reduced mertias about hmge 


' 



Mam body pitch 

slug-ft^ 

641 

631 

623 

Mam body roll 

slug-ft^ 

641 ■ 

631 

623* 

Main body yaw 

slug-ft^ 

3 

3 

4 

Secondary body pitch 

slug-ft^ 

240 

243 

220 

Secondary body roll 

slug-ft^ 

87 

85 

73 

Secondary body yaw 

j slug-ft^ 

327 

328 

293 

Total mertia about hmge 




- 

Pitch 

slug-ft2 

, 881 

874 , 

843 . . 

RoU 

slug-ft^ 

728 

716 

696 ■ 

Yaw ' ' 

slug-ft2 

330 

1 

331 

297- 

Twoaxis hmge 





Pitch sprmg 

ft-lb/rad 

0 182 X lCr2 

0.194 xl(r2 

. 

Roll spring 

ft-lb/rad 

0 382-x 10-3 

0.403 X 10-3 

• 

Pitch damper 

ft-lb-sec/rad 

0.162 

0 155 

• 

Roll damper 

ft-lb-sec/rad 

0.029 

0.0268 


Smgle-axis hmge 


, 



Sprmg 

ft-lb/rad 

i 


0.714 X 10-3 

Damper 

ft-lb-sec/rad 



0.395 

Damper stops 


+27.5“ 

+27.5“ 

+29.5“ 
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since observed variations in their librational behavior may possibly be 
due to differences in their boom properties. The three-axis^ two-body 
gravity gradient stabilization system used on 161 and 163 (ref. 11) 
consisted of three extendable booms arranged in a symmetric pattern about 
the plane of the roll-yaw axes. The primary body was composed of the 
payload and the main boom; the secondary body consisted of the two lateral 
damper booms fixed in a V shape relative to each other. The lateral 
booms were nominally located in the horizontal pitch-roll plane. The 
secondary body was connected to the primary body through a two-axis 
(pitch and roll axes) hinge mechanism employing an eddy current damper and 
a torsion wire spring suspension system. Because of the inherent gyro- 
scopic roll-yaw coupling in the libration of a gravity gradient satellite, 
restriction of the rotational hinge motion of the secondary body to two 
axes is theoretically sufficient to achieve three-axis damping of the 
entire satellite. 

Self-extending SPAR BI-STEM booms manufactured by SPAR Aerospace Products, 
Ltd., of Canada were used on 161 (main boom 1/2" dia., lateral booms 1/4" 
diara.). Self -extending booms manufactured by the Westinghouse Electric 
Corp. of Baltimore, Md. were used on 163 (main and lateral booms 1/2" 
diam.). Both types of booms were interlocked, a feature which tended to 
give these booms a higher torsional stiffness than the open cross section 
booms used on earlier satellites. Perforations were provided on the 
Westinghouse booms in an attempt to better distribute the solar radiation 
energy picked up by the boom and thus reduce the magnitude of thermally 
induced boom distortions. The SPAR booms were not perforated. 

Satellite 164 


The basic geometry of satellite 164 is illustrated in figure 3. The three- 
axis, two-body gravity gradient stabilization system used on 164 (refs. 12 
and 13) consisted of three extendable booms arranged in a symnetric pattern 
about the plane of the pitch-yaw axes. The booms were the interlocked, 
nonperf orated SPAR BI-STEM type used on payload 161. The primary body was 
made up of the payload, main boom and front lateral boom (fixed to the 
payload); the secondary body consisted solely of the lateral damper boom. 
The lateral booms were nominally located in the horizontal pitch-roll 
plane. The secondary body (damper boom) was connected to the primary body 
through a single-axis hinge mechanism that constrained boom motion to a 
vertical plane. The hinge provided hysteresis damping torques and torsion 
wire spring restoring torques. 

The design of this single axis damper configuration was based on the 
inertial coupling concept suggested by Tinling and Merrick (ref. 14). By 
skewing the horizontal principal axis of the secondary body (damper boom) 
out of the orbital plane, all motions become strongly coupled. Under 
these conditions, three-axis damping of the entire satellite is achieved 
by the single degree of freedom motion of the damper boom about its hinge. 
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FLIGHT DATA - ■ 

. 

The three satellites were equipped with attitude instrumentation for 
determining the Euler angle relationships in pitch, roll and yaw between 
the satellite's local vertical coordinate system and its body fixed axes. 

The angles in this case are defined in the usual sense so that in its’ ■ 
preferred equil'ibrium orientation the satellite's body fixed axes are- 
assumed to coincide exactly with its local vertical coordinate system' 

(pitch axis with the orbital angular momentum vector, roll axis withi the 
orbital velocity vector and yaw axis with the ‘local vertical^vector) . 

Attitude Reference System 

The method of solving for these attitude angle’s, described in reference 2, 
depended upon an accurate determination in both local vertical and body 
fixed coordinates of, the direction vectors to the sun and the Earth's . 
magnetic field. ^ A digital computer program, based on tr’adkdng data, was 
used to calculate these vectors in a local vertical coordinate system*; 
satellite sensor data was used to determine these same vectors in a body-, 
fixed coordinate system. A digital computer orthogonal matrix transforma- 
tion was then used to determine the desired Euler angle relationship between 
the two coordinate systems as defined by the two sets of identical vectors. 

The sun data was obtained from a set of three Adcole solar sensors (Adcole 
Digital Solar Aspect System) manufactured by the Adcole Corp. of Waltham, 
Mass. These sensors had a pyramidal field of view of about 128“ and were ' 
judiciously arranged on the top of the payload so that there was nearly 
complete coverage of the celestial sphere. However, certain fields of 
view (e.g., directly over the payload) were not covered, while others were 
covered by two sensors. The sensors measured the angles of the incident 
sunlight with respect to the body -fixed coordinate system of the space- 
craft. These angles, in the form of digital outputs, were sampled and 
stored in the satellite's memory system. 

The magnetic field data was obtained from a triaxial flux-gate magnetom- 
eter (Triaxial Magnetic Aspect Sensor) manufactured by the Schonstedt 
Instrument Co., of Reston, Va. The unit consisted of three sensors orthog- 
onally aligned with the spacecraft's body-fixed axes. Each sensor pro- 
duced an analog output voltage which was dependent upon the magnitude of 
the ambient' magnetic field and the angle between the field vector and the 
sensor's axis. Sampled values of the three output signals, stored in the 
satellite's memory system, were later used to digitally compute the 
direction of the Earth's magnetic field vector in a body-fixed- coordinate 
system. The digital program used in this computation included provisions 
for correcting the flight measurements for that portion of the ambient 


* These calculations were based on an empirical formulation of the com- 
ponents of the Earth's magnetic field obtained by Jensen and Cain (see 
ref. 15) and the known Earth-Sun ecliotic relationship. 
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magnetic field emanating from the spacecraft. This correction was based on 
laboratory pre-flight measurements of the spacecraft's magnetic properties. 

The accuracy of> the attitude reference system depended not only upon the sensor 
resolution and alignment but also upon the accuracy of the tracking data^ the 
computerr- formulation of the Earth's magnetic field and the magnetic field 
compensation. An independent check on the accuracy was obtained by comparing 
the scalar anglei between the sun vector and the magnetic field vector in both 
the local vertical and body-fixed coordinate system. If the error difference 
between these two separate computations was less than 5® the results were 
considered to be acceptable. 

! 

Data Collection 

The flight attitude data provided by NRL covered the first six months of 
satellite operation.- It consisted of printed time histories of each satellite's 
pitch‘d roll and'-'yaw attitudes as determined by the day-to-day interrogation of 
their memory storage systems. The memory systems stored about one day's worth 
of satellite sensor data sampled at a range of about one sample every 154 
seconds*."' These attitude plots were carefully compiled, edited and assembled 
in chronological order. After an initial review it became apparent that 
certain of these plots tended to capture the essential characteristics of each 
satellite's behavior. These individual plots were therefore singled out for 
further examination and have been reproduced in their entirety in Appendix A. 
They have been individually enhanced by tracing through the computed data 
points so as to bring out the distinctive features of each satellite's motions. 
For reference purposes the time for each south-north equatorial crossing, 
starting with Orbit 1 at launch, is also indicated. 



Figure 4 Solar Aspects Data for NR L 160 Series Satellites 


* Real time interrogation was also available for the short period of time 
the satellites were in view of a ground station; during this time, data 
could be directly sampled at rates as high as one sample per second. 
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Solar aspect data for these satellites is shown in figure 4. The percent 
sun in figure 4 is indicative of the period in each orbit in which the 
satellite is not shaded by the Earth's shadow (eclipsed). The sun angle^ 

01 ) in figure 4 is the angle between the sun vector and the normal to the 
satellite's orbital plane. When a=0°, for example, the sun vector is 
perpendicular to one side of the orbit plane, when a=180“ it is perpendicu- 
lar to the other side, and when a=90° it is in the orbit plane. Because ' 
of symmetry, sun angles do not exceed 180°. Since sun angle, a, appeared 
to be such a significant parameter, it has been identified in the 'attitude 
plots in this report by marking the point in time at which a designated 
sun angle was reached . 


SUMMARY OF OBSERVATIONS 

An overall examination of the collected flight data for payloads 161, 163 
and 164 leads to the general observation that there are two basic behavior 
patterns. The satellites are either stable with attitude perturbations 
less than +5° or their behavior tends toward a low-frequency rigid-body 
oscillation dominated by large yaw motions. It is these two' behavior 
patterns that are examined in the following discussion. The events leading 
into and through these patterns are described by referring to figures 5 
to 7. The figures provide information on sun angle versus orbit number 
as well as the location of key events. They begin with the first orbit 
on day 273 of 1969 and end on day 144 of 1970. The broad lines (solid and 
dotted) superimposed on the sun angle line indicate the characteristic 
behavior patterns for those periods when actual flight data for each 
satellite was available. The figures are supplemented in the discussion 
by referring to copies of appropriate sections of the flight data in 
Appendix A and to the ground commands summarized in Table 3. 

Satellite 161 


This satellite displayed a simple form of anomalous behavior. It was either 
very stable or it oscillated in yaw. There were no yaw inversions. ' P,itch 
and roll motions were generally small and they did not appear to play^ a 
significant role in the behavior mechanism. 



♦ Insertion Tronsient 

• Yaw Oscillotions >10* 
o Minimum Point 

— Stable 


Figure 5 Sequence of Events, Payload 161 
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Table 3— Summary of Ground Commands 


Payload 

Orbit 

Command 

Payload 

Orbit 

Command 

161,163 

6 

Primary modulation 

163 

1420 

Thruster 2 off, pitch off 

163 

7 

Boom motor 1 on and off 

161 

1689 

Pitch on 

164 

7 

Prunary modulation 

164 

1689 

PCM bad 

163 

8 

Lateral booms released 

161 

1690 

Pitch off 

161 

20 

Lateral booms released 

161 

1703 

Pilch on 

164 

21 

Lateral booms released 

161 

1704 

Thruster 1 on 

163 

34 

Thrusters 1 & 2 on and 

161 

1709 

Thruster 1 off, pitch off 



off, heaters 1 & 2 on 

164 

1689- 

PCM inoperative 

163 

40 

Heaters 1 & 2 off 


2054 


163 

48 

Command mam boom 

164 

2054 

Last orbit with successful 



m and out 



command 

161,163 

90 

Command memory slow 

163 

2867 

Pitch on 

164 


read 

163 

2868 

Thruster 2 on 

161 

145 

Thrusters and heaters on 

163 

2877 

Thruster 2 off, pitch off, thruster 1 



and off 



on, yaw on due to failure of 

163 

159 

Thrusters and heaters on 



thruster 2 



and off 

163 

2884 

Yaw off 

164 

159 

Pitch momentum wheel 

163 

2889 

Yaw on 



on and off 

163 

2896 

Thruster 1 off 

161 

325 

Voltage control on 





163 

2898 

Yaw off 

161 

335 

Thruster 1 on 

163 

2903 

Yaw on 

161 

394 

Thruster 1 off 

Pitch momentum wheel on 

163 

2908 

Yaw off 

163 

491 

161 

3151 

Heater 2 on to reduce charge 

Thruster 2 on 

161 

494 



current 

163 

494 

Thruster 2 on 

161 

3165 

Heater 1 on for load 

163 

504 

Thruster 2 off 

163 

3165 

Heater 1 on for load 

161 

506 

Thruster 2 off 

163 

3254 

Heater 1 off 

163 

506 

Pitch off 

163 

3262 

Thruster 1 on. 

161 

1096 

Truster 2 on, pitch on 



pitch on 

161 

1099 

Pitch off 

163 

3267 

Thruster 1 off 

164 

1179 

Pitch on 

163 

3268 

Thruster 1 on and off 

164 

1180 

Thruster 2 on 

163 

3270 

Pitch off 

164 

1185 

Thruster 2 off, pitch off 

163 

3276 

Thruster 1 on and off 

163 

1413 

Thruster 2 on, pitch on 

163 

3277 

Thruster 1 on 

163 

1415 

Thruster 2 off 

163 

3282 

Thruster 1 off 

163 

1418 

Thruster 2 on 

163 

3290 

Heater 1 on for load 
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At point see figure 5 and Appendix A^ shortly after insertion into 
orbit and in full sunlight^ the satellite is oscillating as a rigid body 
in response to the insertion transient. The spacecraft has settled down 
into an inverted yaw position. Pitch motion is almost completely damped 
out, roll IS decaying rapidly, and yaw is dieing out slowly. 

Yaw motions continue to fall as the satellite enters into its first period 
of eclipsing orbits, reaching a minimum amplitude of less than +10° at 
around point 2. After this point, however, the oscillations in yaw 
unexpectedly start to grow so that by the time point 3 is reached, the yaw 
oscillations are actually greater than they were at insertion. The yaw 
response frequency during this time was about .73 cycle-per-orbit, varying 
from .75 cycle-per-orbit at insertion to .69 cycle-per-orbit at point 3. 

After the satellite enters its first period of full sunlight (orbit 720) 
the amplitude of this large yaw oscillation starts to decrease and is 
gradually replaced by a low amplitude, decaying one-cycle-per-orbit 
oscillation in yaw, pitch and roll. At point 4 the oscillation has just 
about disappeared and the satellite is extremely stable. After passing 
through the 180° sun angle position (a position in which the vector from 
the sun lined up with the satellite’s pitch axis), the low amplitude one- 
cycle-per-orbit oscillation in yaw, pitch and. roll very slowly reappears. 
This full sunlight region of stability ends as the satellite enters into 
its second period of eclipsing orbits, and the one cycle per orbit motion 
is gradually replaced by the .73 cycle-per-orbit, large amplitude yaw 
oscillation that was seen earlier. 

The oscillation pattern that followed persisted for the next 1300 orbits 
(i.e., until the satellite was again in full sunlight). The behavior 
throughout this long period was repetitive. The yaw oscillation first 
gradually rose in amplitude, then, after reaching a maximum value, it 
slowly fell in amplitude; finally, after reaching a minimum value, it 
started rising all over again. The location of regions of minimum values 
are indicated on figure 5. A good view of this rise and fall pattern can 
be seen by examining the flight data in Appendix A near the minimum value 
at point 5 and the maximum value at point 6. 

Coming out of this long period of eclipsing orbits, the satellite, at 
orbit 2240, enters its second period of full sunlight. As before, the 
yaw oscillations decrease and are gradually replaced by one-cycle-per- 
orbit perturbations in all three attitude traces, point 7. The sun angle 
on this pass, however, does not reach 180°, and the oscillations do not 
completely disappear as they did during the first full sun pass. 

As the satellite leaves full sunlight and enters into a third period of 
eclipsing orbits, the oscillations again begin to grow. The low-amplitude 
one -cycle -per -orb it motions are replaced by the larger amplitude .73 
cycle-per-orbit yaw oscillations, and the rise and fall pattern that was 
seen earlier is repeated. The \Aiole behavior pattern, in fact, is repeated, 
starting all over again at orbit 3160 as the satellite enters into its 
third period of full sunlight. 
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Satellite 163 


The behavior of this satellite was quite complex. It appeared to be sus- 
ceptible to one per orbit pitch oscillations throughout its entire flight 
and large amplitude yaw oscillations and yaw inversions during periods of 
eclipsing orbits. The only prolonged periods of stability were during 

full sunlight. 4 Insertion Transient 

• Yaw Oscillations >10“ 



Figure 6 Sequence of Events, Payload 163 

At point 1^ see figure 6 and Appendix A^ the satellite^ captured in an 
inverted pitch position^ is still responding to the insertion transient. 
Roll and yaw rigid body frequencies are damping out very slowly. Pitchy on 
the other hand^ is responding at a one-cycle-per-orbit frequency and is 
exhibiting no tendency towards dieing out. At point 2 the 180° pitch error 
is successfully corrected by moving the main boom in and out. After this 
controlled inversion, the resulting transient in roll and yaw is still 
damped while pitch continues its sustained one-cycle-per-orbit response. 

As the satellite makes it first entrance into eclipsing orbits, the yaw 
oscillations unexpectedly grow quite large. The satellite rapidly becomes 
unstable in yaw and by the time point 3 is reached, the behavior is so 
erratic that the satellite undergoes a yaw inversion. This undesirable yaw 
performance continues through point 4. In fact, during the entire first 
passage through eclipsing orbits the behavior is marked by numerous yaw 
inversions and several large amplitude oscillations in pitch and roll. The 
yaw frequency in this region seemed to generally be about 1/2 cycle-per- 
orbit . 

The erratic behavior ends soon after the satellite enters its first period 
of full sunlight. By the time point 5 is reached the satellite is very 
stable with no pronounced attitude perturbations. This stability was 
probably maintained throughout the full sunlight period. 
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Although data for the time between orbit 940 and 1450 was lacking^ it seems 
probable that the satellite again became unstable in yaw after entering its 
second period of eclipsing orbits. By the time points 6 and 7 are reached, 
the satellite's yaw behavior is again completely erratic, with several yaw 
inversions occurring near point 7. Large, one-cycle-per-orbit pitch 
oscillations are prevalent, while large yaw oscillations occur that appear 
to be a mixture of one-cycle-per-orbit motions and 1/2 cycle-per-orbit 
motions . 

As the satellite continues into full sunlight the yaw instability once 
again disappears. At point 8 the satellite is again stable, its angular 
perturbations having been reduced to a relatively low level one-cycle-per- 
orbit oscillation in all three attitude traces. In all probability this 
stable characteristic continued until the next eclipsing period was entered. 

Satellite 164 


The behavior of this satellite was markedly different from that of 161 and 
163. It was stable throughout its initial period of eclipsing orbits and 
unstable in yaw during its first excursion into full sunlight. This pattern 
did not persist, however, for during its second passage through eclipsing 
orbits its behavior rapidly deteriorated and the spacecraft ended up in a 
sustained, large amplitude yaw instability. 

♦ Insertion Transient 

• Yaw Oscillations >10® 



Figure 7 Sequence of Events, Payload 164 

At point 1, see figure 7 and Appendix A, shortly after insertion and in 
full sunlight, the satellite has settled down into a well stabilized orbit. 
Its attitude in yaw is inverted* and the perturbations in pitch, roll and 
yaw are small. This satisfactory performance continues through points 2 


* This inversion error was corrected in orbit 159 by energizing and de- 
energizing a pitch momentum wheel. 
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and 3. In fact^ during the entire first passage through eclipsing orbits 
and: first entrance into full sunlight all attitude errors are small. 

During the sequence of events from point 1 to point 3 there were only slight 
changes in attitude behavior; and in general the resulting small amplitude 
perturbations were confined to approximately one-cycle-per«orbit oscilla- 
tLons/in pitch and 1/2 cycle-per-orbit oscillations in yaw. 

As the satellite continues into the 180® sun angle position^ the amplitude 
of the 1/2 cycle-per-orbit oscillation in yaw unexpectedly increases. The 
satellite rapidly becomes unstable and by the time point 4 is reached its 
behavior is completely erratic with numerous yaw inversions and several 
large amplitude oscillations in pitch and roll. Just as suddenly as the 
instability appeared^ however^ it ceases and by orbit 900 the erratic 
behavior has not only disappeared^ but is followed for several days by a 
period of extremely stable operation. 

Entering into the second period of eclipsing orbits the performance of the 
satellite begins to slowly degrade as the amplitude of the 1/2 cycle-per- 
orbit yaw oscillation gradually increases. Shortly before point 5 the 
satellite again breaks into an instability with several successive yaw 
inversions^ ending up at point 5, in an inverted yaw position. What follows 
is a large amplitude^ limit cycle oscillation in yaw at a frequency of 
about 1/2 cycle-per-orbit that persists through orbit 1130 and probably 
longer. 

Data for the time between orbit 1130 and 1280 is lacking. Between orbits 
1179 and 1185^ however^ the satellite was successfully inverted in yaw 
through the use of the pitch momentum wheels so that when the satellite data 
is picked up again at orbit 1280 it is now in a limit cycle yaw oscillation 
about the 0® yaw position. Within the next few days the yaw amplitudes 
become excessive and the satellite's behavior is again completely erratic. 
Yaw inversions now occur nearly every day^ and as typified by the traces 
near point 6^ the instability persists to the very end*. 


EVALUATION OF FLIGHT DATA 

In seeking to provide an insight into those factors that most directly 
influenced the anomalous behavior of the three satellites^ it became 
apparent that any attempt to single out one or two characteristics could 
not easily be substantiated solely on the available flight data. The com- 
plexity of actually defining (at least in a mathematical sense) the inter- 
action phenomenon between, a satellite's internal dynamic properties and its 
external environmental energy sources precludes the simple pin-pointing of 
these critical factors. For example, the effects of aerodynamic torques 
cannot be readily discerned without some additional measurements. Despite 


* On orbit 1689 the PCM telemetry transmission from 164 was lost. Subsequent 
attempts to correct this malfunction or to send commands to the satellite 
were unsuccessful. 
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these restrictions^ however^ the observations sununarized in the previous 
section indicate several relationships that warrant discussion. 

Magnetic Torque 

Although attitude responses due to variations in the Earth's magnetic field 
can be seen in the flight data^ there is no evidence that the generating 
magnetic torques were large enough to contribute adversely to the anomalous 
behavior. For satellite 163 the effect^ although small, seemed to be most 
pronounced in pitch as the satellite passed over the equator, see figure 8. 
Even the magnitudes of these small attitude perturbations are believed to 
be somewhat exaggerated due to computational inaccuracies resulting from 
slight errors in the magnetometer compensation factors and in the empirical 
description of the Earth's magnetic field. 



813 815 820 


ORBIT NUMBER 

ECLIPSE YAW PITCH ROLL 

Figure 8 Magnetic Field Effect on Satellite 163 


An interesting 14-orbit repetitive pattern can be seen in some of the flight 
data (e.g., orbits 606 to 620 on satellite 163). The pattern probably can 
be related to the Earth's magnetic field, since at the end of 14 orbits 
the satellite nearly retraces its path over the Earth's surface. 

Solar Radiation 


The relationship between sun angle and satellite behavior that can be seen 
in the flight data tends to indicate that the effect of either solar pres- 
sure or thermal bending due to solar radiation plays a significant role in 
a gravity gradient satellite's stability. Since a boom's thermal bending 
and twist is related to sun angle, it can be expected that a satellite's 
stability will be influenced by its thermal distortion properties. 

The sudden instability of satellite 164 as it reaches the 180° sun angle 
position (see point 4 on figure 7) is a case in point. The instability 
appears to be related directly with the sun angle, disappearing as soon as 
the sun gets a few degrees away from 180° . 
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Similarly^ in comparing data from 161 and 163^ at least three regions of 
sun-related characteristic responses^ illustrated in figure 9^ can be 
discerned* 

In Region 1 the satellites are in full sunlight with the sun nearly per- 
pendicular to the orbit plane. Both satellites are quite stable with no 
pronounced attitude perturbations. 

In Region 2 they are still in full sunlight but the sun is now inclined 
about 20° to the orbit plane. Although both satellites are still stable^ a 
one-cycle-per-orbit perturbation appears in all three attitude traces. 

In Region 3 the ‘Satellites have gone into an eclipsing orbit with the sun 
inclined about 50° to the orbit plane. Satellite 163 is now unstable in 
yaw with a period of about 2.0 orbits. Pitch perturbations on 163 continue 
at one-cycle-per-orbit. Satellite 161 has begun to undergo relatively large 
yaw oscillations with a period of about 1.5 orbits. Pitch and roll perturba- 
tions on 161 are now small. 

Some perception into the mechanism of thermally induced instabilities of 
gravity gradient satellites can be deduced from Kanning's studies reported 
in reference 6. In this work the behavior of several gravity-oriented 
satellite configurations under the influence of solar radiation was examined. 
The effects of solar pressure torques as well as changing geometry and mass 
distributions due to thermal distortion on the performance of symmetrical 
and asymmetrical satellites were considered for a 1200 km orbit inclined 45° 
to the sun line. Kanning concluded that thermal distortion changes can be 
a critical consideration in the design of an asymmetrical satellite configu- 
ration. Although only one sun angle and only a few isolated configurations 
were examined^ it appears that the simulated performance was significantly 
degraded (none of the cases examined were unstable) by the inclusion of 
thermal distortion. 

An indication of the effect of sun angle on satellite stability can also 
be partially inferred from the recent study by Flanagan and Modi (ref. 16). 
They examined the behavior of a very simple representative satellite (no 
booms) under the influence of solar pressure (thermal distortion neglected) 
and found that in an elliptical orbit the sun angle significantly affected 
the satellite's response. Reviewing their work^ it appears that the 
influence of sun angle would probably be much more pronounced for an 
asymmetrical satellite than for a S3nnmetrical satellite. 

Consideration of thermal bending and twist as a contributing factor to the 
anomalous behavior is complicated to a degree by the effect of "thermal 
twang" (ref. 4).' The "thermal^ twang" excitation is associated with 
eclipsing orbits and is a repetitive disturbance that occurs every orbit. 

In full sunlight the booms are bent due to thermal distortion. Upon enter- 
ing into the earth's shadow the booms return rapidly to an unbent position^ 
introducing an impulse to which the satellite must respond.* The reverse, 
of course, occurs as the booms enter into full sunlight. The shape of the 
resulting impulse can be broken down into a Fourier series so that sustained 
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Satellite 161 


Satellite 163 








satellite inputs can be anticipated at the orbital frequency and its 
harmonics. Since sustained responses during the anomalous performance of 
the three satellites were at frequencies that did not approach these 
"thermal twang" frequencies, it can be assumed that the behavior is not 
being driven solely by these impulses. The impulses on the other hand 
probably contribute to an instability mechanism and may be a feature in a 
feedback path that has not been accounted for. 

Considerable evidence has previously been presented in references 7 to 9 
that relates high frequency flexural oscillations of a sun-lighted boom 
with an instability that is sometimes referred to as thermal flutter. 

There were also some qualitative conjectures in references 2 and 3 that 
the low frequency anomalous behavior of a gravity gradient satellite is 
associated with thermal flutter. Although there is some evidence that the 
final degradation of satellite 164 's behavior was coincident with a higher 
frequency boom oscillation, the preponderance of flight data tends to rule 
out thermal flutter as a controlling factor. 

Dynamic Response Frequencies 

The rigid body frequencies and stability of a hinged two-body satellite 
are controlled to a great extent by the size of the springs used in its 
damper unit. If the springs are too stiff, relative displacement between 
the two bodies is small and energy dissipation due to amplitude dependent 
dan5>ing is negligible. If the springs are not stiff enough, that is below 
some critical value, the satellite will oscillate about a cocked position. 
The ultimate selection usually involves an optimization procedure that 
ends up with springs that have stiffnesses that are slightly above the 
critical value. 

The linearized equations for determining the small amplitude response of 
configurations such as 161 and 163, in which the hinge lies on the principal 
axes, leads to the characteristic equations given in Appendix B. Roll and 
yaw in such a case are decoupled from pitch. The characteristic equations 
for configuration 164, in which the hinge does not lie on a principal axis, 
are given in reference 12 and are somewhat more complex than the equations 
given in Appendix B. Roll and yaw in this case are not decoupled from 
pitch. 

The theoretical response frequencies change as a function of spring stiff- 
ness and in the case of satellite 163 lead to the plots shown in figures 10 
and 11. The plot for satellite 161 is similar. The pitch, roll and yaw 
frequencies for the selected pitch and=* roll springs are appropriately noted 
on these figures. In observing flight data one would expect that the 
frequencies during transient, and even during an anomalous performance, 
would bear some relationship to these characteristic frequencies. This," 
however, was not observed. Instead, as noted on figures 10 and 11, the 
flight frequencies for 161 and 163 (even at small amplitudes)* were close 


* The pitch frequency observed in the flight data for 161 and 163 may be a 
response to a one-per-orbit excitation rather than a transient responses. 
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to the stiff spring frequencies, a condition which cannot occur without 
either assuming some change in the mechanical properties of the satellite 
(e.g., a locked damper spring) or by postulating an additional but unknoxm 
attitude “dependent torque. 




Figure 10. Pitch Frequency Vanation with 
i Spring Constant, Satellite 163 


Figure 11. Roll-Yaw Frequency Variations with 
Spring Constant, Satellite 163 


The observation that the roll frequency on 161 and 163 and the yaw frequency 
on 161 exceeded the theoretical rigid damper frequency is not easily explained 
without assuming some further change in the satellite's structural properties. 
Whether this change was due to an unanticipated variation in inertial 
properties or thermal bending is subject to speculation. 

Because of the inertial coupling in a configuration such as 164, it is 
difficult to ascertain modal responses or to completely distinguish between 
pitch, roll and yaw disturbances. A configuration of this type, because of 
its inherent dependence on coupling for stability, is sensitive to both 
thermal distortion and variation in damper spring properties. Whatever the 
outcome of this coupling mechanism on 164, there is a change in its dynamic 
character, somewhere around orbit 1020, that leads to a rapid degradation 
in its behavior. 
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CONCLUDING REMARKS 


The interaction phenomena obseirved in the flight data appear to establish 
a causal relationship between sun angle and the character of a satellite's 
response. Since thermal distortion and solar pressure are the two 
disturbance factors directly influenced by solar radiation, it can probably 
be assumed that they are influential in the anomalous behavior. Thermal 
distortion properties, bending and twist of each of the satellite's booms, 
may in fact be critical; even though a boom might be perforated to minimize 
thermal bending at normal sun angles, it still may deflect substantially 
at acute sun angles. 

The observations related to discrepancies in the rigid body frequencies 
could be associated with a faulty hinge damper; however, the possibility of 
unanticipated boom deflections inhibiting the operation of the hinge cannot 
be ruled out. Further analytical clarification and classification of the 
mechanism of solar interaction is probably required. If a correct analy- 
tical model of this interaction can be obtained, then the chances ,of 
designing to avoid an instability are much enhanced. 

A case in point is a configuration such as 164 that relies heavily on 
inertial coupling for three-axis stability. Such a configuration is 
particularly sensitive to thermal distortion, and it seems as if such a 
scheme should be avoided until a better understanding of thermal distortion 
effects is obtained. 

A theoretical examination should be directed towards discerning the influence 
of thermal distortion and solar pressure on the long term or orbital 
stability of a satellite as opposed to boom stability. The study should 
consider a complete range of sun angles, boom thermal properties, damper 
unit properties, initial conditions, orbit eccentricity and eclipse times 
consistent with anticipated conditions. 
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Figure A-1. Flight Data, Satellite 161 
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Figure A-1 (continued) Flight Data, Satellite 161 
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Figure A-1 (continued). Flight Data, Satellite 161 
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Figure A*1 (continued). Flight Data, Satellite 161 
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Figure A-2 (continued). Flight Data, Satellite 163 
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Figure A-2 (continued). Flight Data, Satellite 163 
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Figure A-2 (continued). Flight Data, Satellite 163 
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APPEND IK B 


CHARACTERISTIC EQUATIONS 

It can be shown that the small-angle attitude motion of a two-axis hinged^ 
two-body satellite about its equilibrium position is described by the 


following set of linearized equations, 

^]‘^j ^2”^ (B"l) 

n2 + C2’(i?2- Hi) + (12^72- ^2’^! =0 (B-2) 

Cl + Cl” (Cl - C 2 ) + qi«C + “1 Cl - ^1 C 2 = 0 (B-3) 

C 2 ^2 (C 2 ‘ Cl) 92^2? **2 C 2 ~ ^2 (®^) 

? + (l-fj -f2)S2^? -S2f,?j-nf2?2 = 0 (B-5) 


where 

Cj’ = C2/I2, C2’ = C2/I5 

k, = k2/l2. kV = k2/l5 

d, = 3n-(l,-l3)/l2 + k,’ 

d2 = 3n-(l4-I^)/l5 + k2’ 

C,” = Cj/Ij, C2” = C,/l4 

q, = (I,+l3-l2)/Ii 

^2 = (l 4 +' 6 -' 5 )/l 4 

= 4Q“ (1 -q|) + k|/I| 

Ut = -q 2 )+ kj /14 

k] - kj/Ij = ^l/*4 

Hj = pitch motion of primary bodv 
r}^ = pitch motion of secondary bodv 

= roll motion of primary body 
= roll motion of secondary body 


? = yaw motion 

kj = roll spring constant 

k 2 = pitch spring constant 
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C| 5 = roll damping coefficient 
C2 = pitch dampmg coefficient 

- Il , l2» I3 r roll, pitch and yaw mertia of primary body about hmge 
I4, 15, Ig = roli, pitch and yaw mertia of secondary body about hmge 
fi - 2 tt divided by orbital penod 

The pitch equations (B-1) and (B-2) are decoupled from the roll and yaw 
equations (B-3) to (B-5). , The resulting characteristic equation in 
pitch is then 

+ (Cj’ + Cj’) + (dj + d2> s2 + (djC2’ + d2Cj’ - Cj’k2’ - C2’k,’)S 
" • +(djd2-k|’k2’) =0 / (B^) 

The roll and yaw equations (B-3) to (B-5) are all coupled^ a result that 
makes the use of a damper only in roll practical. The characteristic 
equation in roll and yaw is 

+ b4S"‘ + b3S^ + b2$^ + b j S + bg = 0 (B-7) 

where 

bg = -fj -f2) (UjU2 -k2lc2) 

bj = n2(i.fj.f2)(C2”ui+Ci”U2-kiC2”-k2Cj”) 

b2 = -f| -f2)(U| +U2) + UjU2-kjk2 

+ (fjq2lcj + f2qjk2 ^2^2'*!) 

b 3 = n2(,.fj.f^)(d” + C2”) 

+ Cj** (^2 ^^^2^2 " ^2^ 

+C2”(ui + n^f2qi+^^fiqi-ki) 

b4 = .fj + +U2+ n^(f|qj +f2q2> 

bs = Ci” + C2” 
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